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We consider the simulation of interacting high-dimensional systems using pairwise interacting 
qubits. The main tool in this context is the generation of effective many-body interactions, and we 
examine a number of different protocols for obtaining them. These methods include the usage of 
higher-order processes (commutator method), unitary conjugation or graph state encoding, as well 
as teleportation based approaches. We illustrate and compare these methods in detail and analyze 
the time cost for simulation. In the second part of the article, we investigate the influence of noise on 
the simulation process. We concentrate on errors in the interaction Hamiltonians and consider two 
generic noise models, (i) timing errors in pairwise interactions and (ii) noisy pairwise interactions 
described by Master equations of Lindblad form. We analyze and compare the effect of noise for 
the different simulation methods and propose a way to significantly reduce the influence of noise by 
making use of entanglement purification together with a teleportation based protocol. 
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I. INTRODUCTION 

Quantum computers are an example of a wide class 
of quantum-control systems that are capable of simulat- 
ing the Hamiltonian dynamics of any finite-dimensional 
system. In the field of quantum information theory 
such systems are called universal quantum simulators, 
a name first coined by Feynman a quarter of a century 
ago to describe a class of physical systems whose dy- 
namics can be manipulated in order to mimic the dy- 
namics of any other system Feynman suggested 
that such systems might be able to overcome the seem- 
ingly insurmountable problem of efficiently simulating 
quantum mechanical systems. Meanwhile it has been 
shown that universal quantum computers, and all sys- 
tems that can efficiently simulate quantum computers, 
are capable of efficiently simulating the dynamics of all 
finite-dimensional A:-local Hamiltonians and all sparse 
Hamiltonians [!, H[. These classes of Hamiltonians in- 
clude the set of all local spin systems and all Hamiltoni- 
ans that can be efficiently mapped to such systems. 

More specifically, we say that a quantum-control sys- 
tem is a universal quantum simulator on a set of n sub- 
systems if, under ideal conditions, it is capable of gen- 
erating any unitary operation on that set of subsystems. 
There are many physical systems that can satisfy this 
condition given the right forms of control, for instance 
an "always-on" two-body Hamiltonian that entangles a 
set of n qubits can be thought of as a universal quantum 
simulator if is supplemented by arbitrarily fast single- 
qubit unitary control [|[ . These systems can also be con- 
sidered to be universal quantum computers because they 
can efficiently simulate the quantum circuit model (see 
@, 0, 1, S [lfl El El 0, El EE EH for other examples 
involving Hamiltonians manipulated by single-qubit con- 
trol) . This form of quantum simulation is often referred 
to as Hamiltonian simulation, as it involves the manip- 
ulation of a fixed system Hamiltonian. Such universal 
quantum simulators are appealing from an experimental 
perspective because they have the prospect of utilizing 
"global" operations which can sometimes be more easily 
manipulated and created than quantum gates. Other, 
more direct approaches are concerned with simulation 



of specific interaction Hamiltonians (e.g. certain model 
Hamiltonians for high T c superconductivity) using quan- 
tum optical systems such as neutral atoms stored in an 
optical lattice [H EM EH • This form of quantum simula- 
tor, while not universal, are more experimentally feasible 
in the short term. 

Although many universal quantum simulators are 
equivalent to a full scale universal quantum computer, it 
is expected and hoped for that even without using com- 
plex quantum error correction or fault tolerant methods, 
one can simulate specific (relevant) Hamiltonians with 
sufficiently high accuracy, and gain in this way new in- 
sight in the corresponding systems. In particular, quan- 
tum simulators operating on a few tens of spins can be 
expected to be realized on a much shorter timescale than 
a fault tolerant universal quantum computer that needs 
to operate on hundreds of thousand spins. 

With this in mind, in this article we undertake a study 
into the effects of noise on a number of different simula- 
tion protocols. The aims of our investigation are twofold: 

(i) We introduce and investigate different methods to 
generate many-body interactions from two-body 
interactions with help of local control operations. 

(ii) We study the influence of noise in two-body inter- 
action and local control operations on the simula- 
tion process. 

We concentrate on the simulation of many-body inter- 
action Hamiltonians from two-body interaction Hamilto- 
nians. Despite the fact that many-body interactions are 
often neglected in the theoretical description of quan- 
tum systems as they mostly appear as higher order pro- 
cesses, such terms play a pivotal role in the theory of 
quantum simulation. This becomes particularly appar- 
ent when one uses a system of qubits to simulate higher- 
dimensional systems. For instance, if one wanted to sim- 
ulate any eight dimensional spin system with three qubits 
then one would have to be able to simulate three-qubit 
interaction terms. A more extreme example of this is 
shown in [H Hl[ where it is demonstrated that Fermionic 
systems can be simulated using interacting spin sys- 
tems. Such simulations utilize a generalization of the 
Jordan- Wigner transformation [22| to demonstrate how 
Fermionic Hamiltonians can be mapped to a spin Hamil- 
tonian which contains many-body interaction terms. An- 
other application of many-body interaction Hamiltonians 
is in quantum error correcting protocols for adiabatic 
quantum computing. Although it is known that two- 
local Hamiltonians are universal for adiabatic quantum 
computing 23], such schemes are not fault-tolerant. For 
existing error correction protocols it seems that, in or- 
der to correct errors in an adiabatic algorithm that uses 
a two-local Hamiltonian, we are required to use many- 
body interaction Hamiltonians [24[ ■ 

In this paper we particularly focus on the effects of 
noise on protocols for simulating many-body interaction 
terms that are tensor products of the Pauli matrices. 
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Such simulations, while having a simple form, are inter- 
esting because they can often be simulated via protocols 
that can be highly parallelized and they can be used to 
form the building blocks of more sophisticated protocols. 
In addition, such Hamiltonians can transform local noise 
into highly non-local noise. Our approach in this article 
is to examine the effects of noise on the entangling op- 
erations used in these protocols. The noise models that 
we study are quite general and can applied to a variety 
of physical implementations and, in principle, the anal- 
ysis that we perform can be extended to more complex 
simulations. 

More generally, the challenge ahead in the study of 
quantum simulators is to more precisely identify where 
the limitations lie. This problem has two clear pathways, 
the first being to identify the algorithmic constraints on 
quantum simulation. For instance, it is known that simu- 
lating the static properties of quantum systems seems to 
be a computationally difficult task because the problem 
of identifying whether or not a state of a two- local Hamil- 
tonian is the ground state of that Hamiltonian is QMA- 
complete [23j (the quantum analogue of NP-complete). 
Thus it is thought that this task is not efficiently solv- 
able by a quantum computer or any quantum simulator 
that is computationally equivalent to a quantum com- 
puter. In spite of such results, there are currently few 
proven restrictions on the class of systems that are effi- 
ciently solvable on a quantum computer or an equivalent 
quantum simulator. In particular, much work needs to be 
done to identify physical systems that can be efficiently 
simulated by a universal quantum simulator and which 
quantities of these systems can be efficiently extracted. 
For example, recent work [25[ has highlighted that there 
are subtle algorithmic conditions that must be satisfied 
in order to ensure that a quantum simulation protocol 
will achieve any speed-up over a classical simulation. 

The second pathway is to identify the physical restric- 
tions to our capacity to perform quantum simulations, 
which is the main focus of this article. If a quantum 
computer endures too much noise (2r| . its dynamics can 
be classically simulated, the same is true for all quantum 
simulators though the amount of noise that is too much 
for a general quantum simulator is not yet known. One 
way of approaching the problem of how much noise can be 
tolerated by a quantum simulator is to address the prob- 
lem of how to eradicate noise in a simulation. Quantum 
computers can be made fault-tolerant through the use 
of error-correcting codes. The fault tolerance theorem 
tells us that that a noisy quantum computer employing 
quantum error correction can simulate an ideal quantum 
computer without too much additional overhead given 
that amount of noise is below a particular threshold and 
that the errors that occur are of the right kind (for a 
summary see e.g. chapter 10 of [27]). The key message 
of the threshold theorem is that it is in principle possible 
to build quantum simulation devices that can overcome 
the effects of noise and that any fault-tolerant quantum 
computing architecture is an example of such a device. 



Currently, it is not known how far the theory of fault- 
tolerant quantum computing can be extended. The 
first versions of the fault-tolerance theorem applied only 
to architectures that implement the circuit model of 
quantum computing, however, the theorem has since 
been extended to one-way models of quantum comput- 
ing [H, |29|, H3, [U, H2, HH- A s of yet, there is not a 
complete understanding of how to develop a theory of 
fault-tolerance for systems that are controlled by adia- 
batic evolutions and for general time-varying Hamilto- 
nian evolutions. That said, in 34j the authors demon- 
strated that the fault-tolerance theorem can be extended 
to account for systems enduring certain forms of non- 
Markovian noise by considering a model of quantum com- 
putation where circuits are simulated by a "local" time- 
varying Hamiltonian, proving that such Hamiltonian con- 
trol systems are also fault-tolerant. In addition to this 
there have been promising recent developments towards 
developing a theory of quantum error correction for adi- 
abatic models of quantum computing [24| . 

So far we have discussed quantum simulators in the 
context of fault-tolerant quantum computers, yet it is 
far from clear as to whether quantum simulators have 
to be fault-tolerant to be useful. Current fault-tolerance 
thresholds for quantum computers with reasonable noise 
models lie between 10 -5 and 10~ 3 . While theoretical ad- 
vances may raise these thresholds, current belief is it will 
not be possible to build any quantum computer that sat- 
isfies a fault-tolerant threshold for many years. As such, 
it is worthwhile to ask what can be done with systems 
that contain too much noise to be made fault-tolerant and 
yet by many measures are still highly quantum [35]? Such 
systems can simulate classical computations efficiently as 
well as a limited class of quantum systems that one may 
not be able to classically simulate. With the present pa- 
per, where we investigate the influence of noise on the 
simulation of many-body interactions, we aim to shed 
some light on the possibility of using quantum simula- 
tors in such an intermediate regime. 

The paper is organized as follows. In part I, we inves- 
tigate several methods to generate many-body interac- 
tions. In Sec. |TT] we describe the standard commutator 
method that makes use of the Lie- Trotter product ex- 
pansion. In Sec. IIIII we introduce a method that is based 
on unitary conjugation or graph state encoding, while we 
consider teleportation based methods in Sec. IIVI Part 
II of the article is concerned with a detailed study of the 
influence of noise in interactions on the simulation pro- 
cess. In Sec. |V]we discuss distance measures and fidelity 
of noisy processes. In Sec. IVII we describe our first noise 
model, where noisy interactions are described by master 
equations of Lindblad form. We consider a second noise 
model taking random fluctuations in interaction time into 
account in Sec. EH] In Sections EEED El and [X] we inves- 
tigate in detail the influence of noise on the simulation 
process for commutator method, graph state encoding 
and teleportation based method respectively. In Sec. |X] 
we also discuss the usage of entanglement purification to 
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significantly reduce the influence of noise. The results for 
the different methods are compared in Sec. IXI[ and we 
summarize in Sec. IXIII 

PART I: Methods to generate 
many— body interaction 
Hamiltonians 



that is, sequences of applications of the standard Hamil- 
tonian H, together with intermediate local unitary oper- 
ations Ulllk+i, generate (up to higher order corrections 
in St) a unitary operation that is generated by an effec- 
tive Hamiltonian H' (Eq. [1]). Note that more compli- 
cated sequences allow for the simulation of the effective 
Hamiltonian H' with higher accuracy with corrections 
appearing only in higher order St. 



We will consider m systems, each of dimension d, with 
associated Hilbert space H = (C d )® m . In the case of 
d = 2, i.e. qubits, we make use of the Pauli-matrices 
which we denote by <7q = 1, a± = a x ,a 2 = Ci/,C3 = cr z . 
When it is clear from the context, we will often omit 
tensor products, i.e. we identify a\ A ^ <E> = a^<jj B \ 
We will also set % = 1 in the following. 

We assume that the d-leve\ systems interact pairwise, 
and we will be interested in methods to generate effec- 
tive many-body interactions involving up to n < m of 
these systems. For simplicity, we will mainly consider 
d = 2, i.e. qubits. The simulation of arbitrary many- 
body Hamiltonians of m qubits (or, equivalcntly, Hamil- 
tonians of higher dimensional systems with D = 2 m ) can 
be achieved if one is capable of generating 

(i) a specific many-body interaction Hamiltonian H = 
af n for all n < m; 

(ii) fast local unitary control of the individual qubits. 

Hence we will concentrate on the following in methods to 
generate a basic m-body interaction for some fixed m. 

Given (i) and (ii) are fulfilled, standard techniques 
from Hamiltonian simulation can be applied, where inter- 
mediate fast local unitary operations are used to manip- 
ulate the basic Hamiltonian H and generate an arbitrary 
desired effective Hamiltonian 



(1) 



where Hk are m-body interaction Hamiltonians consist- 
ing of Pauli matrices. Here, one makes use of the facts 
that any m-body Hamiltonian can be represented in the 
Pauli basis, and Hamiltonians Hk consisting of Pauli- 
matrices can effectively generated from H via unitary 
conjugation, i.e. 



U k e- UH Ul = e-' ltu » HU K 



(2) 



where local unitary operations Uk suffice. Using the iden- 
tity 



M 



(3) 



we find for sufficiently short times t — St, 



II 



-i8t\ k H k 



II. COMMUTATOR METHOD 
A. Three— body interactions 

A somewhat standard approach to simulating a three 
body interaction from a given two-body interaction is to 
apply a sequence of time evolutions, generated by dif- 
ferent Hamiltonians, each for a short time St. The se- 
quence is chosen in such a way that all first order terms 
in St -when performing a Taylor expansion- cancel, and 
only higher order terms in St remain. These higher order 
terms include products of different two-body Hamiltoni- 
ans, and can hence correspond to effective many-body 
interactions. To be more precise, consider the following 
sequence of time evolutions, generated by the Hamilto- 
nians Hk applied for time St 



j j _ e iH 4 StiH 3 StiH 2 StiH 1 8t 
^tot — ^ ^ ^ 



(5) 



For small St, one can Taylor expand this expression and 
obtains 

4 

U tot = l + i6tJ2 H J 

3 = 1 




(6) 
(7) 



1=2 j=l 

+ 0{St 3 ) 
Considering the case 

Hi = —H 2 , H3 = —Hi, 

one obtains 

U tot = l + ^Cffiffa- i/ 2 i/i)+0(^ 3 ) 

= e ^ 2 x(- 4 /2)[H 1 ,ff 2 ] +0((5t 3 ) _ (8) 

The sequence of interactions described by some properly 
chosen two-body Hamiltonians hence correspond -up to 
higher order corrections- to an evolution which is de- 
scribed by an Hamiltonian that is essentially given by 



-iSt J2k ^kH k 



(4) 



H eB = -i/2[H 1 ,H 2 ], 

which can be an effective three-body interaction. 
Example: For a system of three qubits and 

H 2 = a^aP, 



(9) 



(10) 
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we have that 

e lHuH 2 ]St- +0{st 3 ) ^ 1 (ABC) +i26t 2 a (A) a (B) a p 
+ 0(5t 3 ) 

M e l5 *'^ tf +0(^' 3 / 2 ). (11) 

That is, if we evolve a system as indicated in Eq. ([5]) 
using only two-body interactions with H4 — —Hi, H3 = 
—Hi for a total time of ASt, then the resulting evolutions 
is — up to higher order corrections in St — the same as 
the one resulting from a three-body Hamiltonian 

H eS = a^4 B) 4 C) , (12) 
applied for a time 

St' = 2St 2 . (13) 

Note that there is a dilation factor of 2 /St, i.e. the re- 
quired physical time t to implement an effective three- 
body interaction for a time t' is given by t = 4^/St' /2, 
since the effective three body interaction only appears 
in second order in St. We have assumed that Hx,H2 as 
well as —Hi,— Hz can be implemented. The simulation 
of H c fi is only correct in first order in St' , and unwanted 
terms appear already in order St' 3 / 2 . This corresponds to 
a reduced accuracy as compared to standard Hamiltonian 
simulation schemes for two-body interactions, where the 
desired Hamiltonian is correctly produced up to first or- 
der, and undesired terms (errors) appear only in second 
order. We will refer to this kind of errors as Taylor expan- 
sion errors, which are independent from errors in inter- 
actions and local control operations which will discussed 
in detail in Sec. IIV C II The Taylor expansion errors 
need to be taken into account when using the effective 
m-body Hamiltonian to simulate other Hamiltonians via 
Hamiltonian simulation techniques. 

B. Many— body interactions 

In principle, many-body interactions Hamiltonians for 
arbitrary number of qubits m can be generated in a re- 
cursive way using above method. For instance, one of 
the two Hamiltonians, say Hi, is replaced by an effective 
m — 1 body Hamiltonian. Together with an appropriate 
two-body Hamiltonian a new effective m-body Hamilto- 
nian can be produced. Note, however, that there is a dila- 
tion factor of St/2 in each of these simulation processes. 
That is, the implementation of a m-body Hamiltonian 
for time t' requires a physical time t — 0(t' 2 ( ), if 
basic two-body Hamiltonians are used. 

The time cost can be reduced when using an alter- 
native method, where two n-body interactions are used 
to generate a m = (2n — l)-body interaction. For in- 
stance, the generation of a 5-body interaction has a time 
cost of 0{St 1 / i ), as compared to OiSt 1 ^) when using the 
first method. The desired m-body Hamiltonian appears 
in 0(5t m ^ 1 ) and corrections appear in 0(St m ), where 



m = 5 in our example. Rewriting this in the new effec- 
tive time St m , i.e. one realizes a m-body Hamiltonian 
for time St m , one finds 

St m = 0(St m ~ 1 ), (14) 

and corrections appear in 

0(*C /(m-1) )- (15) 

That is, 

[/• tot = i+tO(*i m - 1 )fl aff + o(ai m ) 

= l + iSt m H eS + 0(St™K m -V) 

w e lStmH ° !! . (16) 

To ensure that the Taylor expansion is a good approxima- 
tion at all instances of the protocol, all involved times, 
in particular St = 0{5t\l <jn ^ 1 ^) (which corresponds to 
the physical time for which two-body interactions are 
applied), need to be sufficiently small. This limits the 
possible values of St m , which is important when consider- 
ing the simulation of general m-body Hamiltonians from 
the basic one. In particular, the total time required to 
simulate a Hamiltonian for time t t ot with Hamiltonian 
simulation techniques (i.e. generating an effective m- 
body Hamiltonian for time St m , and using this Hamilto- 
nian together with intermediate local unitary operations 
to simulate other Hamiltonians for larger times by re- 
peating this process ttot / St m times) requires a total time 
of 

ttot/Stm x 0{5t]l^-^) = t tot O(StJ m - 2 y( m -V). (17) 

The dilation factor (or time cost) 0{Si^ m ~ 2 ^ l ' m ~ 1 ^) ap- 
proaches 0(8t^ L ) for large m and can be significant. Re- 
call that 5tm^ m « 1 needs to be fulfilled, which im- 
plies that the dilation factor, for large m, will typically 
be of order 10 m or larger. 

Example: To make above considerations more con- 
crete, consider the generation of an effective 5-body in- 
teraction Hamiltonian for time St" , 

U tot {St")tte- m " H °«, (18) 

where we use the notation St" = St 5, St' = St%. This 
could take place as follows: (i) use two-body interactions 
UAB^St), Usci^St) generated by two-body Hamiltoni- 
ans ±Hab,±Hbc between systems AB, BC for time 
St = y/St' /2 to produce an effective three-body inter- 
action UABc{St') generated by the effective three-body 
Hamiltonian Habc = ~^/2[Hbc^ Hab]- That is, 

U AB c(St') = U AB (-St)U B c{-St)UA B (St)U B c(St) 

m e Ut ' HABC . (19) 

Use the same method to produce a three-body inter- 
action UcDE(St') generated by the effective three-body 
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Hamiltonian Hcde = [Hde, Hcd]- (h) Use these three- 
body interactions UABc(St'), UcDE(St') to produce an 
effective 5-body interaction U to t(St") generated by the 
effective 5-body Hamiltonian H e g = —i/2[HcDE, Habc] 
for time 

St" = y/W]2= (5t/8)i (20) 

That is, 

U to t(St") = U A Bc(-St')U CDE {-5t')U AB c{St')UcDE(5t) 
w e tSt " H " !{ . (21) 

The choice 

H CD =aP4 D \ U D E = o^af\ (22) 
leads to an effective 5-body Hamiltonian 

H cS = af 5 . (23) 



applied for time t = 7r/4, which is locally equivalent to 
an Ising Hamiltonian H = a z <g> a z . With this choice 
of U, operators of the form ai ® 1 are transformed to 
UpG((Ji(8)l)Up G = ai®cTfu\, with f(i) = if i = 0, 3 and 
/(i) = 3 if i = 1,2. That is, a single-body Hamiltonian 
acting on the first qubit, 

H = Y,*k<7k®t, (27) 

is transformed into a two-body Hamiltonian 

H' = U PG HU PG (28) 

= X a ® (T + AiCTi ® (T 3 + A 2 cr2 <X> CT 3 + A3CT3 <g> fT . 

We remark that the phase-gate corresponds to the uni- 
tary transformation from a product basis \k1k2), ki = 0, 1 
in the a;-basis (i.e. fe« = corresponds to |0) ^ = 
^(|0), + |1),) = ^(|0) + |1)) > and II), = ^(|0)-|1))) 

into a Graph-state basis \G, feife), where G is the graph 
with vertices 1, 2 and a single edge, E 1 = {1, 2}. That is, 

\G,k 1 k 2 ) = U PG \kik 2 ). (29) 



III. UNITARY CONJUGATION AND GRAPH 
STATE ENCODING 

A second method to generate effective many-body in- 
teraction Hamiltonians from basic two-body interaction 
Hamiltonians is by unitary conjugation. That is, before 
[after] the evolution with respect to a single- or two-body 
Hamiltonian H, a (possibly non-local) unitary operation 
U [W] is applied. The resulting evolution is given by 

Ue- itH U* = e - iWHU \ (24) 

i.e. by the transformed Hamiltonian H' — UHlP. If 
U is itself a non-local unitary operation (e.g. gener- 
ated by two-body interactions), then the resulting ef- 
fective Hamiltonian H' can contain many-body terms. 
The unitary conjugation can be viewed as a change of 
basis, where single- and two-body terms of an interac- 
tion Hamiltonian act effectively as many-body terms in 
the new basis. In the following, we will concentrate a spe- 
cific family of basis changes, associated with graph states. 
The corresponding unitary operations U can be efficiently 
implemented using basic two-body interactions, and al- 
low for the systematic construction of the desired many- 
body Hamiltonians via a proper choice of the graph. This 
gives a powerful tool for the construction of many-body 
interaction Hamiltonians. 

Consider as a first example a system of two qubits, and 
a unitary operation U = Upa, 

U PG = diag(l,l,l,-l), (25) 

i.e. a phase gate. A phase gate can e.g. be induced by a 
simple two-body Hamiltonian 

H = {l-a z )®(l-a z ), (26) 



A. Graph— state encoding 

More generally, for systems of N qubits one can define 
an encoding into a graph-state basis corresponding to an 
arbitrary graph G. In this case, the graph G — (V, E), is 
a set V of N vertices connected by edges E, corresponds 
to an interaction pattern (specified by E) between the 
qubits (that are associated with vertices). That is, 

\G,h)= J] UpAk), (30) 

{a,(J}££ 

where k is a binary vector of length N. We define the 
neighborhood N a of a given vertex a as the set of all 
vertices connected to a in the graph, N a — {P\{a,f3} S 
E}. We use the notation 

^=® a&N A a \ (31) 

to refer to an operator acting on all qubits in the neigh- 
borhood of a. Using a such graph state encoding, the 
corresponding unitary operation is given by 

U= II U pg ] - ( 32 ) 

{a,/3}e_E 

Note that the phase gates commute and can hence be 
implemented in parallel. It is now straightforward to 
determine the effect of unitary conjugation with such a 
U. Pauli operators acting on qubit a are transferred as 

follows [ai] 

Ua[ a) U^ = a[ a) a^, 

T T (oOrrt ( Q ) )Vo 

Uai a) rf = 4 a) - (33) 
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We have that a\ a and are transformed into effective 
m-body interaction terms, where m is given by the local 
degree of the graph (i.e. the number of neighbors of a) 

plus one. The action of remains local. The transfor- 
mation of two-body terms follows immediately from Eq. 
([33]), i.e. 

Ua\ a) <g> of V = (Ua^U^iUafu^). (34) 

Note that (Ua^W) and (Uaf ] W) can contain terms 
that act on the same qubit which need to be multiplied. 
Depending on the neighborhood relation of qubits a and 
/?, and on the kind of Pauli operators, the resulting total 
operator can have support on up to |iV Q | + \Np\ + 2 (i.e. 
acts non-trivially on this number of qubits), but may 
also be the identity. For instance, 

= a^a^a^", (35) 

where © denotes the XOR (exclusive or operation), i.e. 
N a © N p = (N a U Np) \ (N a n Np). 

To illustrate the transformation rules, we consider a 
simple example of three qubits arranged as an open chain, 
i.e. with edges E = {{1, 2}, {2, 3}} and correspond- 

( 1 2") f 2 3) 

ing unitary U — U P q UpQ . The two-qubit operator 
transforms to a single qubit operator as 
U<t[ 1} W = a[ 1} a (2) and Uaf } W = af } . Similarly, a two 
qubit operator a[^a^ is transformed to a two-qubit op- 
erator a[^a^ , while a[^a^ transforms to a three-body 

operator cr^cr^ In this case, a three-body opera- 
tor can also be obtained from a single qubit operator, e.g. 
transforms into cr^a^a^ . 

B. Example 1: Three— body interaction 
Hamiltonian exhibiting a quantum phase transition 

We consider now the generation of a three-body in- 
teraction Hamiltonian with transversal magnetic field, 
where we assume a linear chain of qubits with periodic 
boundary conditions. We have 

H = Y,(-< J i a ~ 1) v i i a) 4 a+1) + M a) )- (36) 

a 

This Hamiltonian has been considered in Ref. [3?], [H, H^] 
in the context of optical lattices in a triangular configu- 
ration, where three-body processes may lead in certain 
parameter regimes to an interaction of this form. This 
Hamiltonian exhibits a quantum phase transition. More 
importantly, it has a finite energy gap above its unique 
ground state, a finite classical correlation length but de- 
spite of this an diverging entanglement length [4(| • We 
remark that the ground state of this system for B = is 
a one-dimensional cluster state, i.e. a graph state corre- 
sponding to a closed linear chain. 



We now show how to generate the Hamiltonian H' us- 
ing only two-body interaction and a simple graph state 
encoding. To this aim, we consider a graph with edges 
(2a, 2a + 1). That is, every second qubit is connected 
with its right neighbor. We use this graph for our graph 
state encoding, i.e. define the unitary operation U as 

U = n c/ PG ,2Q+1) - (37) 

a 

It is now straightforward to check the following transfor- 
mation rules 

TT (2a-t) (2a) T rt _ _(2ot-l) (2a) (2a+l) 

Uaf^af^U^ = a{ 2a \ 

Ua {2a) a {2a+1) U^ = a[ 2a+1 \ (38) 
That is, the two-body Hamiltonian 

Hi = -{4^ 1] ^ a) + 4 2a - 1} at ] ) (39) 

a 

transforms into 

H x = 2^~< 7 3 l ^3 ; (4°) 

a 

and gives the desired three-body interaction terms, while 
H 2 = B ^t^?^ + <T? a) °f a+1) ) (41) 

a 

transforms to 

^ = B$>< a) (42) 

a 

and provides the transversal magnetic field in X— 
direction. In total, the two-body Hamiltonian H = 
Hi + Hi is transformed via unitary conjugation with 
U (|3"Tj) into the desired Hamiltonian H (f3"o| . That is, 
for all times t, Ue~ it{Hl+H ^W = e - itH ' . We remark 
that one can in addition achieve a magnetic field in in- 
direction via a single-body Hamiltonian ^2 a a 3 which 
is not changed by unitary conjugation. 

We emphasize that the method proposed here to gener- 
ate the three-body Hamiltonian may be significantly sim- 
pler and easier to implement in optical lattices than the 
original proposal of Ref. [13, [H, HI] , as it does not rely on 
higher-order processes. We require, however, individual 
addressability of the individual qubits in some form. The 
encoding- and decoding pattern (i.e. the operation U) 
can directly be generated from the available two-body 
interaction with is essentially given by a pairwise Ising 
Hamiltonian a z (g> a z resulting e.g. from controlled cold 
collisions. The same interaction, together with fast local 
unitary control operations, can also be used to generate 
the Hamiltonian H = H1 + H2 using Hamiltonian simula- 
tion techniques. That is, the simulation of the evolution 
with respect to the three-body Hamiltonian H' for time 
t to t takes place as follows: 
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1. Use the basic two-body interactions to generate U; 



2. Simulate the evolution with respect to the Hamil- 
tonian H = if i + H 2 for time St by using the basic 
two-body interaction together with fast local uni- 
tary operations. Here, standard Hamiltonian simu- 
lation techniques are applied, and the desired pro- 
cess is approximated up to corrections 0{8t 2 ). 



3. Repeat the process ttot/St times. For St sufficiently 
small the total evolution will be a high-fidelity ap- 
proximation of the evolution generated by H ap- 
plied for time ttot- 



4. Use the basic two-body interactions to generate 
W = U. 



The total evolution after steps 1-4 corresponds to 
e -it t otH ag (j es i re( i. 



C. Example 2: Interacting d— dimensional systems 

We now consider the simulation of Hamiltonians cor- 
responding to (interacting) <i-dimensional systems. To 
illustrate this approach, we first consider a single d- 
dimensional system with d = 2 4 = 16. In contrast to 
the previous example, we will make use of auxiliary sys- 
tems here. That is, we consider a system of five qubits 
1, 2, 3, 4, A and the associated graph G with edges {k, A}, 
k = 1,2, 3, 4, i.e. each of the qubits 1,2,3,4 (which repre- 
sent the 16-dimensional system) is connected to qubit A. 
Here, qubit A serves as auxiliary system and is prepared 
in an eigenstate of a x with eigenvalue +1, i.e. in the state 
|0) x = (|0) + |1))a/>/2. This ensures that any effective 
Hamiltonian that acts on qubit A as either the identity 

or ai A \ leaves this qubit invariant, i.e. qubit A is de- 
coupled from the evolution in this case. However, Hamil- 
tonians involving qubit A can be used to manipulate and 
trigger the effective Hamiltonian acting on system qubits 
1 to 4. The unitary operation U corresponding to the 
graph state encoding is given by 

U=f[U<&\ (43) 
k=i 

This leads to the following transformation rules of 
single- and two-body operators under unitary conjuga- 



tion, where k, ki G {1, 2, 3, 4} 





= *[°\ 


t 


Ua ( 2 a) a (A) tf 


= 4 Q) 


t 




- <7 {a) 

— <7 3 . 


t 


Ua[ ai) a[ a2) tf 


= a[ ai 


L(«2) 

a l 1 




= a[ ai 


1 («2) 

°2 ' 


u4 ai) 4 a2) u^ 


- (7 {ai 

— °2 


L(«j) 
°2 ! 




— °3 


) CT («2) 




(<*2 

= CT 3 






_ _(1)_(2) (3) (4) (A) 
— CTg (Tg (Tg <7 3 O x 



Note that the list contains only operators that leave the 
ancilla qubit A unaltered. As can be seen, all single qubit 
terms, as well as specific two, three and four qubit in- 
teraction terms on the system qubits can be obtained. 
This allows one, in principle, to use Hamiltonian sim- 
ulation techniques to generate any interaction Hamilto- 
nian acting on the four system qubits. Note, however, 
that the generation of effective single qubit unitary oper- 
ations (lines 1-3) requires in part two-body interactions. 
As a Hamiltonian simulation scheme operating on effec- 
tive Hamiltonians may require such effective intermedi- 
ate single-qubit unitaries applied for some time t — 0(1), 
the time cost to simulate certain many-body interaction 
Hamiltonians for time St will be of order t, i.e. a dilation 
factor oft/St. 

In a similar way, two d— level systems of this kind can 
be simulated and coupled pairwise. Consider in addi- 
tion to qubits 1 — 4, A corresponding to a first d = 16 
dimensional system a second set of five qubits 1' — 4', A' 
corresponding to a second d = 16-dimcnsional system. 
We consider for system 1' — 4', A' a similar graph state 
encoding, i.e. the the total graph has additional edges 
(k',A'), k' = l',2',3',4' and U changes accordingly. 
Again, qubit A' serves as auxiliary system and is pre- 
pared in state 10)^, and manipulation of the second d- 
level system works in exactly the same way as for system 
1 (qubits 1-4). In addition, coupling between the two 
d-level systems can e.g. be achieved via the following 
two-body interactions 

UatWW - 4 ai \^A a M A '\ (45) 
U^^Ul = (0 4 =i(7 ( Q ) )(7 K) (7 (^) ; 

The generalization to TV cZ-dimensional systems that in- 
teract pairwise is straightforward. 



9 



D. Example 3: 2D setup with 4— body plaquette 

interaction 

We consider now the simulation of a Hamiltonian that 
contains four-body terms and corresponds to a 2D-setup 
where qubits are arranged on a square lattice. In particu- 
lar, we consider interactions between qubits on the same 
plaquette □, i.e. a Hamiltonian of the form 

h' = Y t 4 a,a) 4 a,a+1) 4 a+1,a) 4 a+1 ' a+1) - (46) 

a 

Such a Hamiltonian can be generated by considering N 
qubits arranged on a rectangular array, and N additional 
auxiliary qubits that are placed in the center of each 
square □ and prepared in the eigenstate of a x , \0) x . The 
graph is such that each of the auxiliary qubits A a is con- 
nected to all system qubits in the corresponding square, 
i.e. to qubits (a, a), (a, a + 1), (a + 1, a), (a + 1, a + 1). 
Hence, single-qubit operations on auxiliary qubit A a 
transforms to a four-qubit operation on the plaquette, 

TT (-Aa)rrt ( Q > Q ) (a-a+l) (a + l.a) (a+l,a + l) A , 

while leaving the auxiliary qubit unchanged. The sim- 
ple Hamiltonian H = B y] a\ Aa \ which may e.g. be 
generated by applying a homogenous magnetic field in 
X -direction, hence transforms to the desired Hamilto- 
nian H ' (|4"6"|) . The construction can be easily generalized 
to higher dimensions (cubic lattice) and other geometries 
(e.g. hexagonal lattices). 

E. Time cost and universal quantum simulation 

1. Fixed graph state encoding: 

The examples considered so far assume that the en- 
coding (or basis change) is fixed throughout the whole 
simulation process. This implies that the time cost to 
simulate a many-body Hamiltonian H = Y2k XkHk for 
time itot , where each of the terms Hk can be obtained by 
a unitary conjugated single- or two-body Hamiltonian is 
given by 

7r/2 + * tot ^A fc . (48) 

k 

Here, we assume that the Hamiltonians Hk can be ob- 
tained directly via unitary conjugation (as in example 
1). That is, one has a constant overhead in simulation 
time, as the basis changes that need to be applied at 
the beginning and the end of the simulation process each 
have time cost 7r/4 (phase gates that can be applied in 
parallel). There is no additional time overhead to simu- 
late many-body Hamiltonians, as the basic single- and 
two-body Hamiltonians are directly transformed to the 
required form. If the method is only used to generate 
a basic m-body Hamiltonian for time St, then the time 



cost is essentially ir/(2St), i.e. essentially the same as in 
the commutator method (Sec. [TTJ) . However, here the 
implementation of the many-body Hamiltonian is exact, 
and the choice of time 5t m , the evolution time of the 
Hamiltonian, is not limited by requirements on the de- 
sired accuracy or the validity of approximations. 

If, as discussed in example 2, the generation of the 
effective many-body Hamiltonian also makes use of 
Hamiltonian simulation techniques -in particular effec- 
tive single-body unitary operations generated by two- 
body interactions-, one encounters an additional increase 
of time cost. Here, however, the graph state encoding re- 
mains fixed throughout the simulation process. 

2. Variable graph state encoding 

Another possibility is to vary the graph state encoding. 
Assume we wish to generate an interaction Hamiltonian 

H = J2*kH k , (49) 

k 

for time short time St, where Hk are interaction Hamilto- 
nians of different kind, e.g. various m-body terms with 
different m. While a particular graph state encoding is 
specially suitable to implement a specific Hamiltonian Hk 
— and in fact such a graph state encoding can be cho- 
sen in a constructive way — , it might be complicated (or 
impossible) to generate all desired terms using the same 
encoding. However, one may change the encoding during 
the simulation process. That is, generate e~ lSt ^ kXkHk 
via a sequence of different unitary conjugations and sim- 
ple single- and two-body Hamiltonians Hk, i.e. 

Y[Uke- iStXk6 "Ul « e - m ^ x » Uk& » u l + 0(St 2 ). (50) 

This is similar to standard Hamiltonian simulation, ex- 
cept that here the intermediate unitary operations Uk are 
not local, but correspond to some graph-state encoding. 

Note that in such a set-up, it is even sufficient to re- 
strict to local Hamiltonians Hk, which are transformed 
via different graph-state encodings (generated by two- 
body interaction) to the desired many-body interaction 
Hamiltonian. That is, we assume that two-body interac- 
tions between any pair of qubits and local control opera- 
tions are available. In this case, the proof of universality 
is very easy. In fact, it is sufficient to show that a proper 
choice of graph state encoding allows one to generate a 
m-body interaction term 

® aeS <t£> ( 51 ) 

on any subset S of parties, where j a € {1,2,3}. Such 
a basic term can be manipulated by means of standard 
Hamiltonian simulation techniques (i.e. intermediate lo- 
cal unitaries, which can be generated from single qubit 
(i.e. local) Hamiltonians and a graph state encoding cor- 
responding to the trivial, empty graph, i.e. Uk = 1) to 
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generate effectively any Hamiltonian acting on the subset 
of S of qubits for some (small) time St. As the subset S is 
arbitrary, the same method allows one to generate arbi- 
trary Hamiltonians, containing terms acting on different 
subsets S. If a Hamiltonian consists of L basic terms of 
the form ([51]) with coefficients Xk, then the time cost to 
simulate such a Hamiltonian is given by Ltt/2 + ^ XkSt. 
A term of the form (p)Tj) is most easily generated by choos- 
ing a graph state encoding corresponding a graph with 
edges (a, a,), where a is some fixed qubit € S, and atj 
runs over all remaining qubits in S (i.e. qubit a is con- 
nected with all other qubits in S). The corresponding 
unitary U transforms single-qubit terms acting on qubit 
a into an interaction term involving all qubits in S, 

which can be transformed via local unitary operations to 
the desired form (|5Tj) . 

The method used in this proof of principle demonstra- 
tion is not 

optimal, as a new graph state encoding is required for 
each individual term. As already demonstrated in two 
examples above, a clever choice of graph state encod- 
ing allows one to generate many terms using the same 
encoding, and hence a significant reduction of the time 
cost can be achieved. In particular, if the Hamiltonian 
corresponds to some translational invariant set-up or ex- 
hibits symmetries, a proper choice of graph state encod- 
ing -where the graph has the same symmetry properties- 
allows in many cases to generate all terms of the same 
form (connected e.g. by the translation operator) simul- 
taneously. For instance, in example 3, two different graph 
state encodings, corresponding to the graph discussed in 
example 3 and an empty graph, are sufficient to generate 
an interaction Hamiltonian with arbitrary 4-qubit inter- 
action terms. This allows, for example, to simulate a ring 
exchange Hamiltonian [4l|, 

H = Y / (b{b2blb 4 + b 1 blb 3 bl)), (53) 
□ 

where the sum runs over all plaquettes □ and opposite 
corners are labelled as 1,3 [2,4] respectively. Note that 
generalizations to other geometries or higher dimensions 
are straightforward. 



to a Hamiltonian of the form 

H = ® aes a ( 3 a \ (54) 

for an arbitrary time t in a deterministic way, where S is 
some set of \S\ qubits. We will consider |5| = m in the 
following. This corresponds to the implementation of a 
unitary operation of the form = exp(—iipH) with H 
given by Eq. [541 and ip = t. 

Different methods to generate specific, non-local uni- 
tary operations are known. Except in special cases (e.g. 
for a CNOT gate), either a sequence of operations using 
a set of entangled states, or entangled states of higher di- 
mension (including some auxiliary particles) are required 
to make the process deterministic. There exist a one to 
one correspondence between arbitrary quantum states E 
(described by a density operator) and completely posi- 
tive maps £ [43j, |4J, |45j . The state E (possibly shared 
between many parties) can be used as a resource to im- 
plement by means of local operations and classical com- 
munication (essentially a teleportation process) the map 
£ on an arbitrary input state probabilistically (i.e. for 
certain measurement outcomes in the teleportation pro- 
cess). In many relevant cases, e.g. for gates of the form 
U v , either or Ul is implemented with probability 
p = 1/2. In this case, the implementation of U v can 
be made deterministic by using a sequence of telepor- 
tation processes [44j corresponding to unitary operations 
U {p ,U2 V , Un v , etc., where additional unitaries are only ap- 
plied if the previous unitary was not implemented suc- 
cessfully. Alternatively, a different entangled state con- 
sisting of additional auxiliary particles can be used for 
teleportation, where the measurement basis of auxiliary 
particles are determined by the measurement outcomes 
of the teleportation process. In 0, |47| such a method 
is demonstrated for the purposes simulating Hamiltoni- 
ans that are diagonal in the computational basis. This 
approach is similar to the measurement based one-way 
quantum computation using cluster states [48| . We will 
discuss both methods in detail in the following. 

Cirac et al [44| demonstrated that the Jamiolkowski 
isomorphism defines a teleportation protocol that allows 
us to use entangled states as a resource for generating 
entangling unitary evolutions. In this section we revisit 
their work and demonstrate several related protocols that 
allow us to use entangled states to generate entangling 
unitary operations. 



IV. TELEPORTATION BASED METHODS 

A third method to generate many-body interaction 
terms is by using a teleportation based approach. That 
is, entangled states are used as a resource to generate by 
means of teleportation an appropriate multi-qubit uni- 
tary operation, corresponding to a time evolution with 
respect to a certain many-body Hamiltonian H . In the 
following we will discuss two teleportation based methods 
that allow one to generate time evolutions with respect 



A. The Jamiolkowski Isomorphism 

We now review some of the features of the 
Jamiolkowski isomorphism. This review is incomplete 
and discusses only those features of the isomorphism that 
we feel are necessary to understand the teleportation pro- 
tocols which will be discussed later in this section. For a 
more complete treatment we refer the reader to [49l [50| . 

The Jamiolkowski isomorphism is a mapping between 
the set of linear operators and the matrix algebra defined 
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on a composite Hilbert space. More precisely, given the 
composite Hilbert space Ti a 1 ®T~Ia and a matrix E defined 
on this composite space then the isomorphism is given by 

(55) 



£(M)=d 2 A tv AA [E A ' A Pi A M A ] 



where H A has the same dimension as H A , M is any 
matrix defined in the matrix algebra of H A (and subse- 
quently is also defined for H A ), and P£ A is a projector 
representing the maximally entangled state 



1$) 



1 



&A 



i) A \i) A 



(56) 



where d A is the dimension of system A, and denote P$ 
|$)($|. The inverse mapping is given by 



E = E / 



(57) 



where X is the identity map. 

It should be clearly noted that the Jamiolkowski iso- 
morphism is defined for the set of all linear operators and 
the set of all matrices for the composite Hilbert space 
TL A ® Ha ■ This set of physical maps is a subset of this 
set of linear operators. Likewise, the set of quantum 
density matrices is a subset of all matrices on this com- 
posite Hilbert space. It turns out, that we can define 
a teleportation-like protocol that helps us to understand 
the circumstances for which the isomorphism is physi- 
cally relevant. 

If we have a maximally entangled state, I*! 1 )" 4 " 4 , and 
we act with a trace-preserving completely positive map, 
£ ®I, on the system then clearly the resulting state, E 
is a quantum state. See figure [1] We shall now see how 
the state E can be used to probabilistically teleport the 
operation £(p) onto some arbitrary state p. Consider the 
composite system A'AA with Hilbert space Ha' ® T~(.A ® 
H.A- Prepare E on the composite system A' A and system 
A is in some arbitrary state, p, 



Pa'AA 



Ei 



(58) 



Then, if we then perform a projective measurement on 
Pa'AA m a maximally entangled basis which includes 
1$)^, the Jamiolkowski isomorphism tells us directly 

that there is a -4- probability that the outcome of this 

A 

measurement will leave system A' in the state £{p). 
Thus, we say that with probability the operation 

A 

£(p) is "teleported to system A'" . Such a process is out- 
lined in figure^ 

1. The Jamiolkowski isomorphism for many-body systems 

We now turn to the case when A is a composite system. 
Following the description in pi9? ] we will consider the 
case where A and A' are many-body systems that have 
Hilbert spaces of the same type. That is we can write, 



n A = n J 



(59) 
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FIG. 1: In order to obtain the state E the completely positive 
map (CPM) £ is applied to system A of the joint system of 
A and A, which is prepared in the maximally entangled state 
Pi A . (Figure taken from Ref. [I~ 



A" 

• ■ 



A,- -A 

^ - _ 

J Projection f~" ^ 



FIG. 2: (Color online) Given the state E on the composite 
system A' and A, the CPM £ is evaluated for an arbitrary 
input state p by taking p as an input at system A. Then 
the joint system AA is measured in a Bell basis containing 

the maximally entangled state P$ A . With probability -k- 

"a 

the desired output state £ (p) is then obtained at system A'. 
(Figure taken from Ref. |4jj) 



where each constituent Hilbert-space has dimension d Ai 



and W 



H A . If we choose the projector P£ A to 



be the tensor product of projectors onto the two-qudit 
maximally entangled states then we can have a 

straightforward way of extending the isomporphism to 
the many-body setting and retain the interpretation in 
terms of a teleportation protocol [50| . More precisely 
is given by, 



pAA 



pAA 



r>A 1 A 1 



(60) 



where 




(61) 



fc=i 



The only difference between teleportation protocols in 
the many-body case and the case that we have already 
discussed is that in this many-body scenario, the tele- 
portation protocol requires N projective measurements 
to be made between the systems A and A instead of one. 
The definition of E in the many-body case is depicted 
in figure [3J The many-body teleportation protocol for 
the implementation of £(p) on system A' demonstrated 
in figured] 
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A A 

FIG. 3: (Color online) In order to obtain the state E the CPM 
£ is applied to the subsystems Ai of the composite system A = 
(Ai, . . . ,An), which are (locally) prepared in the maximally 
entangled states P^ A% . (Figure taken from Ref. [491 ]) 




Input 



FIG. 4: (Color online) Given the state E on the composite 
system A' = (A[, . . .,A' N ) and A = (Ai, . . .,An), the CPM £ 
is evaluated for an arbitrary multipartite input state p by tak- 
ing p as an input at system A = (Ai, . . . , An)- Then the joint 
systems AiAi are (locally) measured in a Bell basis containing 
the maximally entangled state P^ zA \ With probability -h- 

the desired output state £(p) is then obtained at system A' . 
(Figure taken from Ref. [49l ]) 



B. Teleportation-based gates using weakly 
entangled states 

As we have seen, the Jamioikowski isomorphism shows 
us how to create a teleportation protocol that probabilis- 
tically generates quantum operations when given certain 
entangled quantum states. For the purposes of this pa- 
per we are interested in simulation protocols that gen- 
erate unitary operations of the form U(a) — e~ laa * . 
Such a protocol was described in |44| where the au- 



thors demonstrated a non-deterministic simulation pro- 
tocol based on the Jamioikowski isomorphism that suc- 
ceeds with a higher success probability than the proto- 
col that we discussed in the previous subsection. In the 
next section we will demonstrate how this protocol can 
be made deterministic through the use of an additional 
ancillary system. 

1. Quantum control and the Jamioikowski isomorphism 

In order to understand the protocol demonstrated in 
[iij . it helps to have an understanding of the way in 
which control operations on a state, E, are transferred 
to the associated operation £ under the Jamioikowski 
isomorphism. A list of control rules is given in [4^ | , here 
we primarily interested in two of these rules: 

1. Given matrices B\, C\ that can be enacted on sys- 
tem A' and B2, C2 on system A, then the matrix, 

E' = Bf <» B^E £ cf <g> C* 2 A (62) 
is isomorphic to the map, 

£'(M) = B^BlMCDGx (63) 
where the superscript T denotes transposition. 

2. £ is separable with respect to Ak and Ai iff E is 
separable with respect to Ak and Ai. That is the 
CPM corresponding to state E ® F is £ <g> J 7 , where 
£ and T are the maps isomorphic to E and F. 

The first rule can be extended to demonstrate that if 
control is locally limited to parties A and A 1 then the 
allowable set of control operations are limited to local 
operations with one-way classical communication. The 
only allowable classical communication must flow from 
part A to party A! as the manipulations performed by 
party A occur prior to implementing the map £ whereas 
those performed by A occur after. To see this, consider 
the following mixture of matrices, 

E ' = Y1 of R $ E (Qi' ® R t) ] - ( 64 ) 

Using Eq. [H2]we see that the corresponding CPM is, 

£' = Y J Qi£{ R lMR*)Q\. (65) 

3 

If the matrices {Qj} and {Rj} represent measurements 
and unitary operations, then the operation Rj occurs be- 
fore the operation Qj. Thus the set of operations that 
can physically be performed on the state E are restricted 
to, 

E ' = E Qti ® R f E (Qfi ® R f) f > ( 66 ) 
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where the CPM IZp — J2j RjpRj lS bi-stochastic trace- 
preserving CPM 1 , and if j corresponds to a measurement 
outcome at party A the operation Qj(p) = J2i BijpBjj is 
a trace-preserving CPM. It is interesting to note that in 
the case where all operations are either unitary or pro- 
jective measurements then the bi-stochasticity of 1Z is 
assured. This can be seen by noting that for any unitary 
operator U, (U*)^ = U T = (U*)^ 1 and all projectors 
onto pure states are real and Hermitian. 

2. Using weakly entangled states to generate phase gates 

In this section we review the non-deterministic simu- 
lation protocol for generating unitary operations of the 
form U(a) = e~ %aa ^ that was presented in 50]. 

Firstly, we shall give the two-qubit protocol that ap- 
peared in (44j in which it was shown how to generate the 
unitary operator, e - iaa ^ a ^ ) on an y given input state. 
We will then discuss why this protocol works and how to 
extend it to the more general many-qubit case. 

1. Prepare the state, 

\a) = cosa|$ + ) Al ^|$ + ) A3 ^ 

-isma\<5>-) AlA 'i\<5>-) A2A 2, (67) 

where a = ir/2 N , |$+> = 1/V2(|00> + |11)) and 
|$") = 1/a/2(|00) - |11)). The state, |a>, is pre- 
pared between parties A and A' . 

2. Introduce party A which is in the desired input 
state p. 

3. Perform a Bell basis measurements between A\A\ 
and A2A2. Note that we can write the Bell basis 
as 

|* <11 < a > = l®o- il)ia |* + ), (68) 
where ct m = 1, a lj2 = a z , cr 2 ,i = a y , and ct 2 , 2 = 

4. If the measurement on A\Ai produces the outcome 
I ^ii,i 2 ) j a Pply CiiM t° qubit A! 1 . If the measure- 
ment on A2A2 produces the outcome \ ^ , j 1 .j 2 )j apply 
a h,h t° qubit A' 2 . 

5. If i\ = ji then U(a) = e^ 10 " 7 *®"* has been applied 
to system A' and the procedure terminates. This 
occurs with probability 1/2. 

6. If i\ 7^ ji then we have applied U{—a) to system 
A'. This occurs with probability 1/2. In this case, 
we re-label system A' as system A and then prepare 
the state \2a) on A A' and repeat the procedure. 



1 Thus the map TZ = Rj pR*j is also a trace-preserving CPM 



The probability of success, that is that U (a) has 
been applied to system A' , upon repetition is 1/2. 
With probability 1/2 the unitary {/(—3a) is applied 
to system A' . 

7. For the kth iteration we prepare the state |2 fc-1 a). 
After N steps the procedure must necessarily termi- 
nate as U(—(2 N — l)a) — U(a) up to an irrelevant 
global phase. 

Why does this procedure work? It is a variation of 
the teleportation protocol that we saw in the previous 
section. The state \a) can be prepared by the follow- 
ing operation on the pair maximally entangled states 

Looking at the definition of the the Jamiolkowski isomor- 
phism and noting that A, A 1 ', and A are each qubits then 
it is clear that the state \a){a\ = E and it's correspond- 
ing CPM is £ = U{a). According to the teleportation 
protocol of the previous section, in order to apply U(a) 
to an unknown state wc need to perform Bell-basis mea- 
surements on systems A and A. For the state \a) given 
above, a successful run of the protocol coincides with 
measuring systems A\Ai and A2A2 to be in the state 
|$ + ), this occurs with probability 1/16. 

In the protocol that we have presented in this sec- 
tion, we use local unitary operations and one-way clas- 
sical communication to increase this success probability 
to 1/2 for any given instance of the protocol. The key to 
this improved success probability is that the Pauli cor- 
rection operators we apply, er^ j 2 and aj 1 j 2 , allow us to 
manipulate the state a. 

A Bell basis can be generated by applying a single- 
qubit Pauli operator to the state |<1> + ), 

|¥i 1 ,i a > = l®o-i I ,i a |$ + >. (70) 

This equation tells us that performing a measurement in 
this basis and getting the outcome |'I'i 1 ,i 2 ) is is equiva- 
lent to applying the projector |$ + ) followed by the Pauli 
operation <Ji lt i 2 . Interestingly, it is irrelevant which qubit 
the operation <Ji Xl i 2 is applied to. In Step 3 of this pro- 
tocol wc perform Bell measurements between systems A 
and A. If the outcome of these measurements are all |$ + ) 
we have directly applied the gate U(a). Alternatively, we 
can see from Equations (|62p and ([65]) that we will have 
applied the map ^[(cr^ ,i 2 ®Pji,j 2 )p(ci-L,i 2 ® cr ji,h )) ■ ^ n ^ ne 
case where i\ = ji we can apply the corrections of given 
in Step 4 of the protocol above and succeed in apply- 
ing U(a). These corrections work because they commute 
with the unitary U(a). In the case where i\ ^ ji, the 
effect of applying the correction is to apply U(—a) to p, 
as 

U(a)a lul2 <g> a juj2 = a lut2 <8> a jlt j 2 U(-a). (71) 

This protocol works half of the time because U{a) com- 
mutes with half of the set of two-qubit Pauli product 
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operators {aj ® <Jk\j,k = 0, ...,3}. In the next section 
we demonstrate how to make a similar protocol that de- 
terministically applies U(a) to p through the use of an 
ancillary qubit. 

As has been shown in one can generalize this 

protocol to produce evolutions by the n — qubit unitary 
U{a) — e~ laa * . In this case the probability of a success 
in any instance of the protocol is still 1/2 as in the two- 
qubit case. In essence the protocol is the same. Prepare 
the state |a) using n Bell pairs, perform n Bell-basis mea- 
surements, apply n Pauli correction operations. In each 
instance the correction operations are the same as those 
for the two-qubit case. That is, if l^i^a) is measured, 
apply the correction <Ji lt i 2 on the appropriate qubit. Like 
in the two-qubit case, success is in an instance of the pro- 
tocol is achieved in the case when tensor product of all 
of the Pauli correction operations that are applied com- 
mutes with the operator af n , which occurs with a success 
probability of i. 



C. Teleportation-based gates using GHZ-type 
states 

In the previous section we saw how we can use telepor- 
tation as a primitive to stochastically generate unitary 
operations of the form U = e~ laa ® . It has been demon- 
strated [El], H3] that gate-teleportation protocols like this 
are always deterministic if the entangling operation to 
be teleported onto p is an element of the Clifford group 
(defined as the normalizer of the Pauli group) . If an op- 
eration is not in this group, as is the case with U(a), 
then the necessary correction operations in a determinis- 
tic protocol are entangling. 

In this section we shall demonstrate that by perform- 
ing teleportation using GHZ-type states as a primitive we 
can deterministically "teleport" the non-Clifford opera- 
tions like U{a). In essence, by utilizing prior entangle- 
ment with ancillary qubits we can avoid applying entan- 
gling correction operations to the teleportation protocol. 
The simulation protocol that we demonstrate here is es- 
sentially the same as those which was demonstrated in 
the context of graph-state based simulation in [Z(| • 
However, here we use the language of teleportation in or- 
der to highlight the connection with the work presented 
in [Sif. 

Consider the state 

\GHZ) = -L(|+)®«|0> £ + \-r n \l) E ), (72) 

this state is an n+1 qubit GHZ-state. The GHZ state is a 
two-colorable graph state. It is easy to see that this state 
can be generated from phase gates by applying a phase 
gate between qubit E and each of the other n qubits in 
the system. It has recently been shown that such states 
can be purified [H @ H, ^ . 



1. Simulation protocol 

In Section [TV B 21 we saw a teleportation protocol that 
used weakly entangled states as a resource to generate 
U(a) that had a probability of success in any given it- 
eration of 1/2. This protocol failed when the Pauli cor- 
rection operations did not commute with U(a). When 
this occurred, the result of the protocol was to apply the 
operator U(—a) to the state p because the correction 
operation commutes through U(a) to generate U(—a). 
Interestingly, in a failing instance of the protocol, if we 
could "go back in time" and change the weakly entangled 
state | a) to | — a) we would have applied the operator 
U(a) to p. 

In this section we demonstrate an alternate method for 
generating the states \a) or | — a), that is the states corre- 
sponding to the maps U = er laa ^ a - and U = e 1010 *® *. 
This method allows us to generate these states after per- 
forming Bell measurements between systems A and A. 

Consider the following GH Z-like state, 

\k) = ^(|$+) AlA 'i|$+) A2A H0) iS 
v2 

+ l®a z ® 1® cj z \<i> + ) AlA i\<f> + ) MA i\l) E )i73) 

where the superscript E denotes the state of an ancillary 
system. We define the following measurement basis for 
system E: 

{\m),\m±)} = {cos(a)|0) + i sin(a)|l), 

zsin(a)|0) +cos(a)|l)} (74) 
{| — to), I — m±)} = {cos(a)|0) — i sin(a)|l), 

isin(a)|0) - cos(a)|l>}. (75) 

If we re- write the state \k) in terms of the {|to), |to_l)} 
we find, 

\k) = ^[(cos(a)|$+) Al ' 4 'i|$ + ) A2 ^ 
v2 

-i sin(a)l ®<t z ®1® cr 2 |$+) AlA ' 1 \^+) A ^) \ m ) E 

+ (cos(a)l <8> a z ® 1 <g> c7 z \<S> + ) AlA '^ \<5+) A ^ 

+i sin(a)|$ + ) AlA 'i |$+) A ^2) \m±) E ] . (76) 

Which can be simplified to, 

|«) = ^[e- M ^ 1 ' T ^|$ + )' 4lA 'i|$ + ) A2 ^|TO>' E 
v2 

+af 1 af 2 e- la ^' 1 ^' 2 \^ + ) AlA '^+) A2A ^\ni 1 _) E ] 
= ±(\ a )\rn)+afa?*\a)\m ± )). (77) 

Eq. [77] demonstrates that a measurement on | k) of sys- 
tem E in the {|to), |toj_)} basis will either generate the 

state | a), if the outcome is to, or the state at 1 a A ^\a) , 
given toj_ as the outcome. By performing a similar anal- 
ysis it is possible to see that a measurement in the basis 
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{| — to), | — mi)} produces either the state | — a) or 

the state a A± a^ 2 1 — a) . Given that we can always per- 
form Pauli correction operations on system A' after this 
measurement, we see that a measurement on the ancilla 
system E can generate the states \a) or | — a). 

With this technique in mind, we shall now outline our 
new protocol. 

We first consider the two-qubit case. We want to 
deterministically generate an evolution by the unitary 
U = e - %aa *® a z on the unknown two-qubit state p for 
any a. 

1. Prepare the resource state \k) 

2. Introduce party A which is in the desired input 
state p. 

3. Perform Bell basis measurements between A x Ax 
and A 2 A 2 . 

4. If the measurement on A\A\ produces the outcome 
l^ii.ia), apply (Ji lt i 2 to qubit A\. If the measure- 
ment on A 2 A 2 produces the outcome J2 ), apply 
<Jj 1 j 2 to qubit A' 2 . 

5. If i\ — ji, then measure system E in the 
{|m), |toj_)} basis. If the outcome is m the proce- 
dure terminates. If the outcome is m± then apply 
a <r z operation to qubits A\ and A 2 to terminate 
the procedure. 

6. If i\ 7^ ji, then measure system E in the {| — m), | — 
m±)} basis. If the outcome is to the procedure 
terminates. If the outcome is — m± then apply a 
a z operation to qubits A[ and A' 2 to terminate the 
procedure. 

7. Upon termination, system A' will have undergone 
the evolution 

£{p A ') = U(a)pU( a y. (78) 

Essentially, the procedure that we have outlined here 
is the same as that for weakly entangled states. The key 
difference is that in the case of the protocol for weakly 
entangled states, we constructed |a) first. Under that 
protocol, we are restricted to manipulations that are lo- 
cal operations with one-way classical communication. In 
this protocol, we keep systems A and A' entangled with 
system E until the final step. The effect of measuring 
system E is to collapse the total composite system into 
the desired state, either \a) or | — a). 

This protocol is easily extended to generate evolutions 
of the form U = e~ laa ? n . Firstly prepare the state 

\k) = ^(|$+) AlA 'i...|$+) A " A "|0) 
v2 

+af\..af"\<£ + ) A i A K..\<t>+) A ^\l)). (79) 



Then perform the same Bell measurements and Pauli cor- 
rections as above. If the operations applied after the mea- 

A' A' 

surement in Step 4 commute with the operator a z 1 ...a z " , 
then measure in the {|m), |m_i_)} basis and apply the op- 

A' A' 

erator a z 1 ...a z " if the outcome is m±. If they do not, 

measure in the {|— to), |— mj_)} basis and apply at 1 ...af" 
if the outcome is toj_. 



PART II: Influence of noise 

As we discussed in the Introduction, building a fault- 
tolerant quantum computer is a challenging task that 
may not be accomplished in the near future. In the 
meantime it is hoped that non-trivial quantum simula- 
tion experiments could still be performed. In particular, 
the demonstration of an imperfect quantum simulation 
that could still outperform the best-known simulation al- 
gorithms on a classical computer would be a significant 
experimental achievement. 

In Part II of this paper, we discuss the influence of 
noise and imperfections on the simulation of many-body 
Hamiltonians using the different methods described in 
Part I. We focus our attention on protocols for simulat- 
ing Hamiltonian evolutions of the form H a = e~ lacr ® m . 
As we saw in Part I, while Hamiltonians such as H a have 
a particularly simple form, they can be used alongside 
the Lie- Trotter formula to simulate the dynamics of any 
finite-dimensional Hamiltonian. Our focus is on these 
more simple Hamiltonians because in many cases the ef- 
fects of noise on these protocols can be analytically stud- 
ied. This gives us insight into the types of errors that 
might occur in more sophisticated simulation protocols. 
In particular, the following analysis could be used to de- 
rive error bounds on the measurement statistics of more 
complex simulations that use Hamiltonians like H a as a 
building block. 

We begin Part II by discussing general methods for 
classifying noise processes. We then go on to introduce 
the noise models that we consider, local dephasing and 
depolarizing noise and errors induced by random errors 
in the evolution time of entangling Hamiltonians. We 
then set about integrating these noise models into the 
protocols that we discussed in Part I and examining the 
effects analytically, where possible, and computationally. 
Mostly, we look at the effect of introducing noise to the 
entangling operations in these protocols and consider lo- 
cal operations to be perfect. This is done because for the 
purposes of quantum simulation we are concerned with 
the way that the noise is spread throughout the system 
as a result of the entangling operations in each protocol. 
In addition, we often focus on regimes where the effect 
of local noise on single-qubit operations can be incorpo- 
rated into the noise model for the entangling operations 
without much loss of generality. 
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V. DISTANCE MEASURES AND FIDELITY 

We will use two different measures to quantify the qual- 
ity of the simulation process. Time evolution with re- 
spect to a specific many-body Hamiltonian H for a time 
t corresponds to a unitary operation Ut (which is the 
desired evolution), while time evolution governed by a 
master equation of Lindblad form (describing the noisy 
evolution), leads to a completely positive map £ t . We 
consider the case where the operation is acting on m 
qubits and determine how close £t is to Ut- 

We use two devices for comparing the closeness of £ t to 
Ut, Jamiolkowski fidelity [58j and the local noise equiv- 
alent. The mathematical detail of both of these devices 
are discussed in more detail below. It is worth noting 
here that both the Jamiolkowski fidelity and the local 
noise equivalent have useful operational interpretations 
when being used to describe a simulation protocol. 

A. Jamiolkowski fidelity 

The Jamiolkowski fidelity F(Ut, £t) measures the close- 
ness of one quantum operation to another. This fidelity is 
based on the Jamiolkowski isomorphism that was intro- 
duced in Section IIV Al that relates quantum operations 
to quantum states. The Jamiolkowski fidelity measures 
the overlap between two quantum operations by mea- 
suring the overlap between their corresponding states as 
defined by the Jamiolkowski isomorphism. This will be 
more rigorously defined below. 

There are many ways that one might wish to go 
about defining fidelities for comparing quantum opera- 
tions. The Jamiolkowski isomorphism is useful because 
it is easy to calculate, it satisfies some important math- 
ematical properties, and it has several clear operational 
interpretations which are closely linked to the theory of 
quantum computation and simulation [58| . 

The first operational interpretation of the 
Jamiolkowski fidelity is in terms of a sampling quantum 
computation. Imagine that we were to implement a 
quantum operation, U, on some initial state |a;) in order 
to attain measurement outcomes, y, on the output state. 
Say that we vary the input x and sample in order to 
attain the joint probability distribution p x (y). Such a 
scenario is the way that one might perform any number 
of physically relevant simulation protocols. Say that 
the ideal distribution, that which occurs if U is imple- 
mented perfectly, is p x (y) and the actual experimentally 
achieved distribution is q x (y) , which occurs as a result of 
implementing and imperfect operation, £. Then it can 
be shown [58J] that the Jamiolkowski fidelity provides a 
lower bound on the overlap between the two probability 
distributions. That is, 

F(U,£)<B(p x (y),q x (y)) (80) 



stances are chosen uniformly at random. The Bhat- 
tacharya overlap is a measure for comparing probability 
distributions, while it doesn't have a succinct operational 
interpretation, it can be used to bound the Kolmogorov 
distance between two probability distributions. 

Let us now consider a different simulation protocol, 
one that is based on a function computation. Imagine 
that we desire to use a unitary process, U , to simulate 
an evolution from the state \x) to the state \ f{x)). We 
then perform measurements on the output state \f(x)) 
in order to determine the property of that state given 
by f(x). If we instead performed the evolution £ and 
output the state £(\x}(x\), then the average probability 
p e , of getting an incorrect value of f(x) given a uniform 
distribution over x is bounded by, 

p e <l-F(U,£). (81) 

More succinctly, if an experiment is designed to simu- 
late some system in order to deterministically evaluate 
some property of a state, then the Jamiolkowski fidelity 
provides a bound on the average probability of an error. 

We now move on to the more mathematically oriented 
properties of the Jamiolkowski fidelity. There are a num- 
ber of ways of defining a fidelity for a quantum operation 
that are based on the Jamiolkowski isomorphism. In the 
case where we are always comparing an ideal unitary op- 
eration to a CP map it is convenient to use the following 
definition: 

F{£ u U t ) = F(E t , |tf t » - (* t \Et\Vt), (82) 

where 

i* t ) = iw®ui B) m, 

E t = £ t (B) |$) AB ($|, (83) 

and |<f>) = ^-/V2™Y^k=Q X W)A\k)B is a maximally en- 
tangled state of 2m qubits, where we labelled the basis 
of m qubits by {|0), |1), . . . |2 m - 1)}, i.e. identified the 
m qubits with a 2 m -dimensional system. The operation 
acts on one part of the maximally entangled state. This 
fidelity was called the Jamiolkowski process fidelity in 
[58j ] , as it can be only requires elements of an operation's 
process matrix in order to be calculated. In this paper 
we simply refer to it as the Jamiolkowski fidelity. 

It is worth noting that the Jamiolkowski fidelity is lin- 
early related to the average fidelity, F(U,£), between a 
unitary operation U and a CP map £ , 

F(U,£) = J d^\U^£(\i>)(i;\m), (84) 

where the integral is over the pure states \ip). For a sys- 
tem of dimension d, the average fidelity can be calculated 
from the Jamiolkowski fidelity via, 



where B(p x (y), q x {y)) is the Bhattacharya overlap be- 
tween the probability distributions and the problem in- 



(85) 
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B. Simple bounds on the Jamiolkowski fidelity for 
simulation protocols 

The Jamiolkowski fidelity, as defined for this paper, 
can be used to generate several distance measures on the 
the set of all quantum operations. One such measure, 
D(£t,Ut), will be used in this paper [58|, 

D(S t , U t ) = y/l-F(£ t ,Ut). (86) 

This distance measure is a metric on the space of quan- 
tum operations that is stable (that is remains unchanged) 
under the addition of ancillary quantum systems. Impor- 
tantly, this measure also satisfies the following chaining 
inequality, 

D(U 2 o U x ,£ 2 o£{) < D(U u £x) + D{U 2 ,£ 2 ), (87) 

where U± t 2 are unitary operations and £\ t2 are trace- 
preserving quantum maps. This chaining inequality ba- 
sically tells us that if we compare two sequences of quan- 
tum operations, then the distance between the two se- 
quences will not diverge more rapidly than the distance 
between each individual step. This inequality turns out 
to be useful in our context, as quantum simulation pro- 
tocols can often be cast in terms of sequences of unitary 
operations. 

The chaining inequality allows us to obtain an up- 
per bound on the distance (and correspondingly a lower 
bound on the fidelity) between two control sequences. 
Given an ideal unitary sequence U = U m o f7 m _x ' ' '°Ui 
and the corresponding noisy sequence £ — £ m o £ m -\ o 
• • • o £ i then the chaining inequality tells us, 

m 

D(U,£) <Y,D{U 3 ,£ 3 ) < mD max , (88) 
i=i 

where D max = maxj D(Uj,£j). 

Consider now the case where we wish to compare an 
ideal operation U with a noisy operation of the form, 

£ = Uo£ u = £ u oU. (89) 

where £jj commutes with U . Using the definition of 
F(U,£) it is easy to see that, 

F(U,£) = F(I,£u). (90) 

Correspondingly, we find 

D(U,£) = D(I,£u). (91) 

While this seems quite obvious, such circumstances occur 
readily in physical situations. For instance a system with 
Hamiltonian proportional to o~f m undergoing dephasing, 
timing errors (as we will see later), or alternatively for U 
and £ that are both generated by fast control pulses (as is 
the case when using the Lie- Trotter summation formula) . 

Each of the simulation protocols discussed in this pa- 
per use sequences of unitary operations followed by mea- 
surements. The ideal unitary operations that appear 



in each protocol are either local unitary operations (of- 
ten referred to as LU operations) or are applications of 
the entangling Hamiltonian a z ® a z between pairs of the 
qubits in this system for some time t, that is the unitary 
Ut. In the case where we can consider all LU operations 
to be perfect it is very easy to establish upper (lower) 
bounds on the Jamiolkowski distance (fidelity) using the 
the chaining inequality. 

A generic control sequence that appears in this paper 
looks something like, 

U = (LU) m o U t o (LU) m -! o U t o ■ • ■ o (LU) 2 o U t o (LU) 1: 

(92) 

where (LU)j denotes some local unitary operations and 
Ut is an fixed entangling unitary that could in princi- 
ple be applied between any set of qubits in the system. 
If we consider noise that only occurs as a result of the 
entangling operation Ut then we can replace Ut in the 
above sequence by some quantum map £ t that describes 
the noisy entangling operation and our noisy control se- 
quence looks like, 

£ = (LU) m o£ t o (LU) m -i o £ t o • ■ ■ o (LU) 2 o£ t o (LU)i. 

(93) 

Comparing these two sequences with the chaining in- 
equality we find, 

D(U,£) < D{(LU) m ,{LU) m ) + D{Ut,£t) 

+D((LU) m -!, (LU) m -i) + D(U t ,£ t ) + ■■■■ 

(94) 

If we note that for all j D((LU)j, (LU)j) = then the 
above expression can be simplified to, 

D(U,£)<kD(U u £t), (95) 

for a sequence involving k implementations of the map 
£t ■ Thus the problem of establishing an upper bound on 
the distance between U and £ is reduced to finding the 
distance between U t and £ t . The corresponding fidelity 
bound is: 

F{U,£)>l-k 2 D{U u £tf (96) 

This bound is not the tightest bound that can always 
be established analytically. Later in this paper we shall 
see similar constructions that are optimized depending 
on on the simulation model that appears and the noise 
that is being modelled. Better bounds can be found by 
grouping the basic operations in each simulation control 
sequence that commute with each other before applying 
the chaining inequality. For instance, if all noise terms 
and unitary operations in a control sequence commute 
with one another the fidelity can be calculated precisely. 
Conversely, if each local unitary operation that appears 
does not commute with Ut then it is not possible to find a 
better bound on the fidelity without resorting to further 
approximations. 

While the Jamiolkowski fidelity is in general a suitable 
distance measure for noisy quantum operations, one has 
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to be a bit careful about its interpretation. To be precise, 
when considering unitary operations U t with small t, i.e. 
time evolution with respect to the desired Hamiltonian 
for a small time t, the operation will be very close to the 
identity. The influence of noise, as resulting from a mas- 
ter equation description we consider here, also increases 
with time t. For small t, the influence of noise is hence 
very small, and also the noisy operation will be close to 
the identity. This will lead to a fidelity close to 1, also 
in cases where the noise dominates over the Hamiltonian 
part of the evolution. In particular, it can happen that 
for small t the fidelity is close to one, but the noisy op- 
eration is not capable of creating entanglement (i.e. the 
corresponding operator Et is separable [44j]). 



C. Gate fidelity for long— time evolutions 

To allow for a better comparison, we thus also con- 
sider unitary gates generated by time evolutions for long 
times of O(l). The ideal evolution is given by the uni- 
tary gate Ut, which is obtained by the time evolution 
with respect to the desired Hamiltonian H for time t, 
where e.g. t = tt/4. The ideal evolution is compared 
with a sequence of basic noisy evolutions that approxi- 
mate the desired gate. To be precise, the noisy evolution 
Est for (possibly small) times St is applied n times, where 
ri = t/St, to approximate the desired unitary operation 
Ut- The Jamiolkowski fidelity 



F t (£ St ) = F(Y[£ Su U t ), 



(97) 



fe=i 



is a suitable quantity to measure the quality of Est and 
compare such noisy evolutions with different times St. 
Also from a practical point of view, F t is a useful quan- 
tity, as in many cases the goal is to generate a certain 
unitary gate, and F t measures directly with which fidelity 
this is possible. It turns out that the optimal strategy is 
often not to generate very precise approximation of the 
desired evolutions for short times St and apply a sequence 
of many of such evolutions, but to chose some fixed inter- 
mediate time St' which only require a moderate repetition 
of these basic evolutions. 



VI. NOISE MODEL 1: MASTER EQUATION 

We will describe the influence of noise by a master 
equation of Lindblad form that is derived under certain 
assumptions, e.g. the Markov approximation. To be pre- 
cise, we model the interaction between systems by a mas- 
ter equation of the form 



9 <w 



Cp 



(98) 



where we have separated the evolution into two parts, a 
Hamiltonian part, 



Tip :— —i[H, p] 



corresponding to the ideal unitary process in question, 
and a noise part described by some Liouvillian superop- 
erator C. The formal solution of this master equation 
can be written as 

p{t) = e^ +£ >V(0), (100) 
which can, for small t = St, be Taylor expanded 



p(St) 



P(0) 

st(H + r)p(o) 
st 2 



(H +C + CH + H£)p{0) 



0(St 3 



(101) 



In the following, we will often assume a coupling of the 
individual particles to independent thermal reservoirs, 
that is noise acts independently on particles. The Li- 
ouvillian C is in this case given by a sum of terms corre- 
sponding to different particles, 



(102) 



For qubits, we describe the coupling of a single particle 
to a thermal reservoir by a quantum-optical noise term 

M, 



D 



r>(a) 13 / 1 \r (a) (a) (■>) (mi . ini m; 

L> > p = — — {l — s)[al r 'a y _'p + pa + 'a_'-2a y _'pa + '\ 



r («) fa) 



" r ai oi . (a) lay r> (a i 

——s[a_ 'cr + ' p + pa_ 'a K + ' — 2a\ 'pa\ 



2C-B 



[2p-2a^pa<?\ 



(103) 



with a^ = l/2(ax ±io~y^) act on particle a and 2C > 
B. While parameters B,C give the decay rate of inver- 
sion and polarization, s £ [0, 1] depends on the temper- 
ature T of the bath. More precisely s = lim^oo ( 1 + cr ~ ) t , 
where s = 1/2 corresponds to T = oo. For B = and 
C = 7, this corresponds to dephasing noise, while for 
s = 1/2 and B — C = k, we obtain white noise, corre- 
sponding to a depolarizing channel. Without interaction, 
i.e. for H — 0, the corresponding time dependent com- 
pletely positive maps obtained by integrating the master 
equation are given by a dephasing channel, 

V^p = p(t)p+ 1 —^(p + a^pa^) , (104) 
with p(t) = e (t or a depolarizing channel, 



M {a) p = p(t)p + 



1 ~P(t) \ - fa) fa) 



E^rvr, (105) 



3=0 



(99) with p(t) — e Kt respectively. 
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A. Local noise equivalent 

We also use a third quantity to determine the quality 
of the simulation process, which is the local noise equiv- 
alent (LNE). Here, the noisy process described by S t is 
compared with a (virtual) noisy process S t that corre- 
sponds to an evolution with respect to a superoperator 
given by 

m 

n + J2 c(a) ( 106 ) 

a=l 

for time t, where the Hamiltonain part Tip = —i[H,p\ 
corresponds to the ideal evolution, while the noise part 
is described by individual (local) coupling of qubits to 
thermal baths (see Eq. I103[) . In particular, we consider 
the cases of dephasing noise (B = 0, C — 7) and white 
noise (s = 1/2, B = C = k). The parameter 7 [k] 
determines the strength of coupling of each qubit to the 
reservoir, and is adjusted in such a way that the processes 
St and St — exp[(H + J^^Li £( a ))i] lead to the same 
Jamiolkowski fidelity with respect to the ideal unitary 
evolution Ut, i.e. 

F(S t ,U t ) = F(£ t ,Ut) (107) 

The coupling strength 7 [k] for which Eq. 11071 is 
fulfilled hence give the local noise equivalent, i.e. the 
amount of local noise that would lead to same fidelity 
after an evolution for time t. This provides a suitable 
measure for the quality of a noisy evolution. We remark 
that typically St and St will not be identical. Neverthe- 
less, if the initial process containing two-body interac- 
tions is described by single qubit dephasing/depolarizing 
noise with some parameter 70 [kq] , a LNE 7 > 70 [k > kq] 
can be considered as increase of noise level, where it is 
implicitly assumed by such a comparison that generation 
of a many-body interaction is as difficult as the gener- 
ation of a two-body interaction (which might, however, 
not be the case). As we will see in the following, some 
methods typically lead to a LNE which is much larger 
than the initial 70, Ko- 

In certain cases, the LNE can be determined analyt- 
ically from the Jamiolkowski fidelity. Consider e.g. the 
case where the ideal (desired) evolution is given by a 
Hamiltonian H = erf™, and we consider dephasing noise 
on individual qubits described by a Liouvillian C^ a 'p — 
—j/2(p — a z a ^ ' puz^)- In this case, the Jamiolkowski 
fidelity of the process S t — exp[(H + ^ a ^ a ')t] ~ 
exp(7it) exp(^ Q £( a H) is given by 



F{£ t ,Ut) = F((g> a v( a \l) = 




where 2?( Q ) is the dephasing channel defined in Eq. 11041 
with p(t) — exp(— 74), and we used that Hamiltonian 
and noise parts commute. The relation of the LNE to 



the fidelity F = F(St,Ut) of a process is in this case 
given by 

ln(2F 1 /" 1 - 1) , s 

7= K - '-. (109) 

For other Hamiltonians, and depolarizing noise, where 
Hamiltonian and noise parts do not commute, it is often 
difficult to solve the corresponding master equation an- 
alytically. However, the LNE can be easily determined 
numerically also in these cases. 



VII. NOISE MODEL 2: TIMING ERRORS 

As we have discussed, all of the Hamiltonian simula- 
tion methods in this article can be thought of as being 
implemented via control sequences of Hamiltonian evo- 
lutions. In these protocols we manipulate the amount of 
time that a Hamiltonian evolution is induced on a sys- 
tem. In this section we examine the effect of inaccuracies 
in this evolution time. 

Say we wish to physically implement the unitary oper- 
ation U = e~ lHt and we give ourselves control over the 
amount of time that the evolution occurs. In reality, we 
can only control t with finite precision. In a single in- 
stance of a controlled evolution by H for some time t, we 
assume that there will be some Gaussian distributed ran- 
dom error 5 with mean 0. The corresponding completely 
positive map is given by, 

/DO 
d5p(S)e- iH ^pe iH ^+ 5 \ (110) 
-00 

where p(S) is a Gaussian distribution. 

Importantly, an error such as that described by Eq. 
1 1 1 01 can be thought of as two distinct, commuting op- 
erations on the the state p. That is, we can define 
two super-operators U t (p) = e~ lHt pe lHt and T{p) = 
JZo d5p(5)e~ ms pe lHS = dSp(S)U s (p) such that, 

S(p) = ToH(p)=HoT(p). (Ill) 

Noting this, we are free to think of the map given in 
Eq. 11101 as being comprised of some "perfect" unitary 
operation, that is the operation U, which is followed by 
a noisy operation, T. 



1. Unitary Hamiltonians 

An evolution according to Eq. 11101 has a particularly 
simple form in the special case where a Hamiltonian is 
of the form H = ojH u , where = 1. This property 
allows us to express the unitary evolution generated by 
such Hamiltonians in the following simple form, 

e - lHt = cos(wt)l - i smiwtfHu. (112) 
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There are a number of interesting Hamiltonians in this 
class. For instance the two-qubit Ising interaction, as 
well as single-qubit Pauli rotations. 

Utilizing the identity in Eq. 11121 we can express the 
map T(p) as, 



T(p) = / d5p(5) [cos 2 \wS)lpl + sin 2 (cj5)H u pH u 

J — OO 

-i sin(wJ) cos(ujS)(H u pl - lpH u )] . (113) 

As p(5) is an even function of S and the product 
sin(w<5) cos(luS) is an odd function, we know that, 



p(S) sm(u8) cos(w<5) = 0. (114) 

' — OO 

Thus wc find, 

/OO 
dSp(S) cos 2 (oj5)lpl + sin 2 {uj8)H u pH u . 
-OO 

As we have assumed that the timing error is Gaussian 
distributed around zero we can write the probability dis- 
tribution p{$) in the following way 



P(S) = 



1 



CTV27T 



(116) 



where a is the standard deviation of the probability dis- 
tribution. We can find an analytic expression for T{p) 
by inserting Eq. 11161 into Eq. 11151 using some standard 
trigonometric identities, and noting the following integral 
identity, 



(117) 



dxe ax cosibx) = \ l — e 
u 4a 



Doing this we find that the completely positive map asso- 
ciated with a random Gaussian error in the control time 
of a Hamiltonian is given by, 

T(p) = l(q + l)lpl+±(l-q)H u pH ul (118) 



where q(a, to) 
as, 



-la u 



Equation [TT~8l can be re-written 



T{p) = qp+^-Jl{p + H uP H u ). (119) 

which has a similar form to that previously seen for the 
de-phasing channel. 

The corresponding Jamiolkowski fidelity a for random 
timing error on a unitary Hamiltonian would be, 



F(H,£) = F{1,T) 



(1+9) 



The corresponding Jamiolkowski distance is 
D(H,£) -- 



(1-9) 



(120) 



(121) 



Example: Consider the two-qubit Ising Hamiltonian, 
H = uja z (£> a z . If wc were to evolve a two-qubit system 
by H to a system of two qubits for some time t, then 
the map associated with a random timing error on the 
evolution time would take the form: 



T(p)=qlpl+V—ti(p + a z 



a z pa z ®a z ). (122) 



Such a map can be thought of as a collective dephasing 
channel. We shall use this example frequently in the 
remainder of this paper. 



2. Sums of commuting unitary Hamiltonians 

We can extend the results of the previous section to 
a wider variety of Hamiltonians. Consider the following 
Hamiltonian, 



H 



(123) 



where each Hj is a unitary Hamiltonian, uij is real, and 
any two Hj commute. As each of the Hj commute, an 
evolution according to such a Hamiltonian can be factor- 
ized into a sequence of n unitary evolutions, 



or in terms of super-operators, 

u t {p)=i4o...oi4(p). 



(124) 
(125) 



The effect of a timing error on such an evolution can 
be analyzed in a fashion similar unitary Hamiltonians 
that we examined in the previous subsection. We can 
describe an evolution of H with a timing error using the 
completely positive map £ that was defined in Eq. IllOl in 
the previous section. Like we saw in the previous section, 
we find that because all of the Hj commute £ = Ut°T~H = 
Tjj o U t . Hence in order to understand the effects of the 
timing error we only need to consider the nature of the 
operator Th . 

In terms of superoperators, Th can be written, 



Th{p) 



p{5)U}o...oUUp), (126) 



where the order the sequence of U( is irrelevant. Note 
that the integral in the expression for Th is only over a 
single variable, 8. As we did in the previous section, we 
can expand each of the U 3 S terms in Th using Eq. 11121 
and integrating out odd terms in the expansion. That is, 



MP) 



p{5)U l s 



(cos 2 (w n J)lpl + sin 2 (uj n 8)H n pH n ).(l27) 
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Some algebra demonstrates that if each U J S is expanded 
in this way, then the resulting expression becomes quite 
complicated and writing it here would be unhelpful. Gen- 
erally, the expression for Tjj involves a sum of all prod- 
ucts of the unitary Hamiltonians Hj . The nature of such 
an expression depends greatly on what each of the Hamil- 
tonians Hj actually are. This is best illustrated by an 
example. 

Example: Imagine we wish to apply an Ising Hamil- 
tonian to 3 qubits of the form: 

H = uJa l z a 2 z + (128) 

for a time t. The resulting timing error operation can be 
written, 

/oo 
d5p(5)(e" ~ /m ' 

-OO 

This can be expanded in the fashion described above, 

/oo 
cWp(£)[cos 4 (w5)]lpl 
-OO 

+ cos 2 (w(5) sin 2 (ujS)(a l z a 2 z pala 2 z + ala 3 z pala 3 z ) 
+ S m i (cjS)a 2 z a 3 pa 2 z a 3 z ]. (130) 

This expression can be simplified to, 

T H (p) = gtpt + (1 - g)(o\(j 2 z po\o 2 z + a 2 z a 3 z pa 2 z a 3 z 

+ (1 - g falalpalal (131) 



(129) 



where 



d5p{5) cos 4 (w5). 



(132) 



Assuming that p(6) is a Gaussian distribution with mean 
it can be written as it is in Equation 11161 and we find 
an analytic expression for g, 

.9=| + V W + ^( e - W ) 4 - (133) 

In this example, the Jamiolkowski fidelity is simply, 
F(H,£)=g. 

From this example we see for any H that is a sum of 
commuting unitary Hamiltonians that it is possible to 
derive a lower bound on the Jamiolkowski fidelity given, 
a random timing error. Expanding T~h we find that there 
is always a term with the following form: 

/oo n 
dS P {S)Y[cos 2 {u Jj 5)l 0n pl 0n . (134) 
j=i 

This means that the Jamiolkowski fidelty must be at 
least, 

/oo " 
d5p(S)Y[cos 2n {cj J S). (135) 

3=1 



If each ujj — lj and p(S) is a Gaussian, then it is possible 
to find an analytic form for this lower bound by using 
basic trigonometric identities, 



F(H,£) > 



1 (In 

22" 



1 



22n-l 

k=0 



(136) 

As our three-qubit Ising model example demonstrates, 
this lower bound on the fidelity can be saturated. Satu- 
ration of this inequality occurs when no product of the 
Hamiltonians Hj are equal to the identity. That is, the 
only term in 7# which gives a superoperator proportional 
to 1 is Ih ■ There are a many Hamiltonians that have this 
property. Two examples include 1-D Ising chains with 
open boundary conditions and Ising interaction Hamil- 
tonians whose interaction pattern forms a star graph. 

The Hamiltonian H = ■_ 1 uJjHj can be simulated by 
a sequence of n independent evolutions of the Hamilto- 
nians Hj . If each of these evolutions had a timing error 
then Eq. 11261 would have n independent Sj terms over 
which one would integrate. This is equivalent to apply- 
ing n different error maps like that given in Eq. 11191 
That is the timing error map would be given by a prod- 
uct of maps like Th 1 ••■'?ir„ ■ The fidelity of a timing error 
in such a case is then lower bounded by 



F(H,£) = F(l,T Hl ...T Hn 



3 = 1 



1 



(137) 



This bound on the fidelity is saturated for the same H 
as the bound derived for Tu above. 

If we consider the situation where the above fidelity 
bounds are saturated, we can compare the effect of tim- 
ing errors on an evolution according H alone versus the 
effect of simulating H = $^?=i wHj by n evolutions ac- 
cording to its constituent Hamiltonians Hj . For simplic- 
ity, we consider the case where u>j = tu and the standard 
deviation in the timing error for implementing H is the 
same as for implementing each of it's constituent terms 
Hj. The resulting fidelities are given by the following 
integrals; 



F(1,T H ) = 



d6p(S) cos 2n (ujS) 



(138) 



for the case where we implement H directly; and 

(poo \ n 

J dSp(5) cos 2 (uS) J (139) 

for the case where H is implemented by evolutions ac- 
cording to the Hamiltonians Hj. We have reverted to 
writing these expressions in their integral form because 
in this form it is simpler to see how to apply Holder's 
inequality to them. Recall Holder's inequality, 



\C(x)ri(x)\dx < / \C(x)\ r dx 




r,(x)\ w dx 



(140) 
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where C( x ), v( x ) are integrable complex functions and 
r, w > 1 satisfying 



1 1 



= 1. 



(141) 



We set ((8) = pfijV, ,,(,$) = p (5)i cos 2 (luS), r = ^ 
and w = n and substitute these values into Holder's in- 
equality 2 . Then, noting that p{5) and cos 2 (toS) are both 
always positive we find that 

)H />oo 
< / d5p((5)cos 2Tl (^). 
— GO 

(142) 

Thus we find that for 5 = Sj and LOj = lo for all j 

F(1,T H ) > F(l,T Hl ...T H J. (143) 

VIII. COMMUTATOR METHOD 

We will now analyze the influence of noise, described 
by some master equation of Lindblad form, on the sim- 
ulation of many-body interaction Hamiltonians gener- 
ated by using the commutator method. The commutator 
method is based on a sequential application of evolutions 
for times St with respect to different Hamiltonians in such 
a way that, when performing a Taylor expansion in St, 
all first order terms vanish and the second order terms - 
that include the commutator of the initial Hamiltonians 
and hence many-body interaction terms- remain. The 
analysis is usually done at the level of operators (see Sec. 
Hi)) , where one obtains (up to higher order corrections) 
an effective Hamiltonian given by H c s — — i/2[Hi, H2] 
which appears in second order St. One can formally per- 
form the analysis also at the level of superoperators, as 
is required when considering noisy interactions described 
by master equations. 



A. Three— body interactions 

We consider a sequence of four operations, each applied 
for time St, where we denote the Hamiltonian part of 
the j operation by Hj , and the Liouvillian part by Cj . 
An initial state p(0) has evolved under this sequence of 
operations after a total time t = 4St as follows 



Pit) 



_ (Hi+Ci)St {H 3 +C a )8t (n 2 +C 2 )St (Hi+Ci)8t 



p(0). 



(144) 



2 M.B. would like to thank Otfried Giihne for pointing out this 
choice of variables. 



For small St, Eq. 11441 can be Taylor expanded using 
Eq. 11011 and one obtains 



P(t) 



0=1 i=i 



St 



4 1-1 



, 3=1 1=2 3=1 



St 



T E E 2 ^ L i + L{H i + AA)P(O) 



1=2 3=1 



St 



2 

0{8t s ) 



3=1 

3\ 



(145) 



As for ideal evolutions, we consider the case where 
Hi = -H a , H 3 = -Hi, i.e. H 4 = -H 2 ,H 3 = -H x . 
In this case, Ylj=i ^jP(O) = 0' an d tne second line of 
Eq. pS] simplifies to 2St 2 Hp{0) with Hp = -i[H,p] 
and H = ^[Hi,H2\. Also some of the terms in lines 
3 and 4 of Eq. 11451 are simplified. Note, however, that 
the term StJ^^i ^j'P(O) -which is first order in St— re- 
mains. While Hj correspond to coherent processes that 
can annihilate each other, Cj correspond to incoherent 
noise processes, where no such interference is possible. 
This implies that noise appears in first order St, while 
the (desired) Hamiltonian part only appears in second 
order, St 2 . We denote additional noise terms of lines 3 
and 4 that appear in second order, St 2 , in Eq. 11451 by 

and use the short-hand notation L = J2j. 
The total time evolution can hence be written as 



r . 
-1 3 ' 



Pit) 




^)+0(St 3 ) 

P (0) + 0{St 3 ) 

p(0) + 0(St 13 / 2 ). 

(146) 



where we have introduced, as in the noiseless case, a di- 
lated time St' — 2St 2 . That is, when applying the se- 
quence of evolutions specified in Eq. 1 1441 for a total time 
of ASt, this is equivalent (up to higher order corrections), 
to an evolution with respect to a Liouvillian superopera- 
tor 

r ri 

(147) 



for a time St' — 2St 2 . Note that the Hamiltonian part of 
the evolution is given by H = | [Hx , H 2 ] which provides 
the desired many-body interaction term, but the noise 
part of the evolution, dominated by £/V2St' is enhance 
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by a factor of {25t')^ 1 / 2 , or equivalently (2St)~ 1 , as com- 
pared to the two-body evolutions. This is a significant 
increase of noise level, which clearly limits the applicabil- 
ity of the commutator method. In particular, the increase 
of noise level becomes larger for smaller times St, while 
on the other hand smaller St increase the precision with 
which the desired many-body interaction terms are gen- 
erated. Recall that the commutator method is based on 
validity of Taylor expansion, which requires small St and 
precision increases with decreasing St. In other words, 
the dilation factor in simulation time given by the time 
cost, translates in the case of noisy interactions directly 
into an increase of noise level. 

Example: We consider in the following a simple ex- 
ample that illustrates the increase of noise level. We 
consider the case where, independent of the interaction 
Hamiltonian, the Liouvillian describing the noise part 
of the evolution corresponds to single-qubit dephasing 
noise, i.e. is given by 



££(«V = £ 7o /2(p-4<VW), 



(148) 



for all particles a that are involved in an interaction. 
For Hamiltonians Hi , H2 given in Eq. [10] and H3 = 
— Hi,H.4 = —H2, the effective evolution resulting from 
the sequence of evolutions specified in Eq. 11441 corre- 
sponds -in the ideal case- to an evolution with respect 
to the effective three-body Hamiltonian 



H c ff = a 



- 



(149) 



for time St' = 2St 2 . When Hi is applied the correspond- 
ing Liouvillian is given by d = L (A) + C (B) . If St < 1 
and 70 << 1 then the evolution Hi + Ci corresponds 
to an evolution according to the Hamiltonian Hi and a 
single-qubit dephasing channels acting on qubits A and 
B. The resulting fidelity is given by 



Fi 



l+p 



(150) 



where p — e~ loSt . Assuming the same fidelities for each 
of the two-qubit operations used to simulate H c g we can 
use the simple fidelity bound derived in the previous sec- 
tion, Eq. [96] to give a lower bound on the fidelity of the 
simulation 



F > 1 - I670&. 



(151) 



Note that this bound is only valid for small St and 70 . In 
terms of the effective simulation time St' it reads, 



16 

F>l-- 70 > 



"Sf 



(152) 



We improve on this bound below by directly calculating 
the fidelity under similar approximations. 

The effective Liouvillian C e s describing the noise part 
of this simulation is dominated by the term —fi= = 

J V2St' 



^%=T, where C x = C 3 = + & B \C 2 



£(B) + £(C) as on i y qu bits AB or BC are involved in 
a particular operation. We thus obtain 

2L^ + 2£( c ) + 4£( B ) 



(153) 



In principle, the sequence of operations leading to the ef- 
fective Hamiltonian H c g can be symmetrized by permut- 
ing the roles of qubits A, B, C, either probabilistically 
or sequentially, which leads to symmetrized noise in the 
effective Liouvillian 



£ c ff, 



8/3(^)+£( B ) +£( c )) 



sym 



(154) 



where the factor 8/3 results from (2 x 2 + 1 x 4)/3, i.e. 
in two out of three cases we obtain a factor 2 in front of 
the Liouvillian term (role of particles A,C), while in 1 
out of 3 cases there is a factor 4 (role of particle B) - see 
Eq. 11531 The LNE with respect to dephasing noise can 
be directly read off from above equation and is (approx- 
imately) given by 

(155) 



where the expressions are only valid for St' -C 1. The 
Jamiolkowski fidelity can be estimated via Eq. 11081 and 
one obtains 

fl , -v / ^ 7 8 7 o/(3v / 2)\ 3 , 

F « ( «1- 2V2Sf l0 . (156) 

We remark that dephasing noise is in many cases not an 
appropriate model, and also the assumption that dephas- 
ing noise occurs independent of the Hamiltonian might 
often not be fulfilled. Consider e.g. the case where a 
basic two-body Hamiltonian H = <j z a z is manipulated 
by means of local unitary operations to simulate a new 
effective two-body Hamiltonian, e.g. H' — <J x a x . In this 
case, the dephasing noise is also transformed, and one 
obtains bit-flip noise for the new effective Hamiltonian 
instead. 

A similar calculation can be performed for depolarizing 
(white) noise, i.e. the Liouvillian describing the noise 
part is a sum of local terms, ^ Q £S a \ with given by 
Eq. 11031 with s — 1/2 and B = C = n. Again, one finds 



- T Ko, (157) 



3V2St' 

for the LNE, while for the fidelity we can simply bound 
it in the limit of small k and St using Eq. [96] by 
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F > 1 



V2 



(158) 



and we can explicitly estimate the fidelity in the same 
limits to be 



F 



1 + 3 e -V6FK, /{&y/2) \ ' 



1 - 3\/»k . (159) 
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Above calculations should be considered as estimates 
of LNE and fidelity. A complete calculation involves also 
additional errors due to Taylor expansion and noise terms 
occurring in order St' . We have thus performed numer- 
ical simulations taking all errors into account, and con- 
sidered the (often) more realistic model of white noise 
(depolarizing channels). The results of these simulations 
for the generation of an effective three-body interaction 
are shown in Figs. 15161 Plots correspond to coupling 
strength k that would lead to single qubit white noise 
channels with parameter p — 0.9; 0.99; 0.999 respectively 
when applied for time St = tv (see Eq. 11051) . The results 
fully agree with the analytic estimates (see Eqs. I157I159[) 
for small St'. For larger 5t' , higher order corrections due 
to Taylor expansion also play a role and further increase 
the noise level (and hence increase the LNE and reduce 
the fidelity). We remark that although for small St' the 
fidelity is close to unity, the quality of the simulation 
is nevertheless very low as the local noise equivalent is 
enhanced significantly. 
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FIG. 6: (color online) Increase of Local noise equivalent of 
three-qubit interaction exp(— iSt'af 3 ) applied for time St' 
generated by commutator method. Plot shows k/ko, as a 
function of St'. Curves from bottom to top correspond to 
p = 0.9 (blue), p = 0.99 (red), p = 0.999 (green), where 
p = exp(— 7tk). Dotted line is the analytic estimate given in 
Eq.H5ZI«/«o = i7 ij ? . 

tion IVffl (Eq. 1 1 20[) we explicitly calculated the resulting 
fidelity of a timing error for such Hamiltonians. Consid- 
ering such a noise model and given that St « 1 then 
we can use Eq. [SB] to estimate a lower bound on the re- 
sulting fidelity. Assuming additionally that the standard 
deviation, a, of each timing error is the same we find the 
following lower bound on the fidelity of the simulation of 
H c ff to be 

F> 1-8(1-?) (160) 

where q — e~ 2a . 



FIG. 5: (color online) Error of three-qubit interaction 
exp(— iSt'af 3 ) applied for time St' generated by commuta- 
tor method. Plot shows deviation from ideal evolution mea- 
sured by the fidelity, 1 — F, as a function of St' . Curves from 
top to bottom correspond to p — 0.9 (blue), p — 0.99 (red), 
p — 0.999 (green), where p = exp(— ttk). Dotted lines are 
analytic estimates of Eq. 11591 1 — F = 3V2St' «o- 



B. Timing errors 

In Section fVIII we derived CP maps that describe the 
effect of random timing errors in the implementation 
of a Hamiltonian evolution. As we saw in Section [TT] 
we can simulate H c g — a z ® a z ® a z by a control se- 
quence involving the Hamiltonians ±ifi = ±a z A ^a z B ^ 
and ±H 2 = ±<Jy <r z . Both H 1 and H 2 are unitary 
and Hermitian, thus they are their own inverses. In Sec- 



C. Many— body interactions 

It is straightforward to perform a similar analysis 
for many-body interactions involving more than three 
qubits. Essentially, one finds that the time cost of sim- 
ulating many-body Hamiltonians using basic two-body 
interactions translates to an increase in noise level. As 
shown in Sec. IIIB1 the desired m-body Hamiltonian ap- 
pears in St m = OiSt^™ 1-1 )), where we have denoted by 
5t m the time for which the effective m-body Hamiltonian 
is applied. The noise, however, still appears in first order 

St, or equivalently 0(Stm'" m This corresponds to an 
increase of noise level by a factor of order 0(<5t~( m ~ 2 )), 

or equivalently 0(Stm m ~ 2 m )■ The total evolution 
for time St m is thus essentially governed by a Liouvillian 
superoperator of the form 

H + 0{5t7 m {m - 2)/{m ~ 1) )C, (161) 
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where 7i = — i[H e {i, p] and H c h is the desired m-body 
Hamiltonian, while C is given by an appropriate sum 
of Liouvillians describing noise in two-body interactions 
involved in the process. This noise is enhanced by a factor 
of 0(8t m {m ~ 2),{m - 1] ). 

The LNE also increases accordingly, i.e. 

7 «O(^ m ( ™- 2 )/(™- 1 )) 70 , (162) 

while the fidelity is approximately given by 

F « 1 - m 7o x 0(8t-< m - 2 M m -V). (163) 

In above expressions, we have again assumed a simple 
model with only dephasing noise specified by a coupling 
constant 70 [7] respectively. For large m, the dilation fac- 
tor approaches 8t~^ and can be significant. In particular, 
the resulting CPM may no longer be capable of creating 
entanglement for all times. 

We have also performed a full numerical analysis tak- 
ing all errors and imperfections into account. We find 
that qualitative scaling for small simulation time 5t m 
is essentially as described above. For larger St m addi- 
tional errors due to Taylor expansion and/or noise terms 
in higher order St further increase the noise level. 



D. Simulation of strongly entangling multi— qubit 
gate 

Above considerations are concerned with the genera- 
tion of a basic m-body Hamiltonian of the form H c g — 
af m for time 5t m . However, such a basic m-body Hamil- 
tonian only serves as a building block to generate - 
via standard Hamiltonian simulation techniques- other 
many-body Hamiltonians, or multi-qubit gates. To il- 
lustrate the accuracy of such simulation in the presence 
of noise, we consider two examples: 

(i) The simulation of the three-body Hamiltonian 



Hi 



_ „(1)„(2)„(3) 



<Jz &z for time t = 7r/4; 



(ii) The simulation of the three-body Hamiltonian 



H 9 



(1) (2) (3) 

cri cri 'cri 



El 



&x f° r time t = 7r/4; 



Example (i) only contains the basic Hamiltonian, while 
example (ii) also includes single-qubit terms and the ap- 
plication of Hamiltonian simulation techniques to gener- 
ate sequentially the two non-commuting terms af 3 and 
Efc=i ^ • Again, we determine fidelity and LNE as a 
function of the time 5t m for which the basic two-body 
interaction Hamiltonian is generated. Small times 5t m 
require several application (and additional intermediate 
local unitary operations) of the basic 3-body Hamilto- 
nian to achieve an evolution for time t = 7r/4, while 
larger times St m require fewer applications but involve 
larger Taylor expansion errors in both the generation of 
the basic 3-body Hamiltonian and the simulation of the 
desired Hamiltonian H from H a s . The results of the nu- 
merical simulation for the Hamiltonian (ii) are shown in 



Fig. [7] for different values of the white noise parame- 
ter p, where p — exp(— itk). That is, the corresponding 
coupling strength to the thermal bath n is such that the 
evolution with respect to the noise part of the Liouvillian 



for time t = n leads to a depolarizing map (see Eq. I105[) 
with parameter p. One observes that there exists an op- 
timal time 5t m that leads to a maximal gate fidelity, and 
the optimal value of 8t m depends on p. Recall that the 
LNE to generate the basic 3-body hamiltonian increases 
with decreasing 8t m (see Fig. |5J), while St m cannot be 
chosen too big to limit Taylor expansion errors. The op- 
timal value of 8t m is hence a compromise between the 
two competing requirements. The results for simulation 
of the Hamiltonian (i) are very similar and are not shown 
here. 
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FIG. 7: (color online) Error of three-qubit gate 



fc=i 



ffc) 

cri gener- 



exp(-ZTr/4#) with H 2 = a\ L> arar + £ 
ated by sequential applications of exp(— i8t m o^ o^) and 
exp(— iSt m 53t =1 <y^). The three-qubit interaction is gener- 
ated by the commutator method, and the figure shows the 
derivation from ideal gate as measured by the fidelity, 1 — F, 
plotted against different values of 8t m . Curves from top 
to bottom correspond to p — 0.9 (blue), p — 0.99 (red), 
p — 0.999 (green), where p = exp(— tik). Notice that there is 
an optimal choice for 8t m to maximize the gate fidelity (see 
main text for details). 



E. Simulation of quDits using qudits 

We have illustrated the commutator method for qubits, 
however the results can without further effort be directly 
adopted to quDits. When using pairwise interacting d- 
level systems to simulate an arbitrary Hamiltonian of a 
single D-level system with D > d, an effective m-body 
interaction Hamiltonian is required, where 



logD 
logd 



(164) 
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and \x\ denotes the next largest integer. This imme- 
diately follows from the fact that the Hilbert space of 
m interacting qudits has dimension d m , and D > d m 
is required that the m qudits can simulate a single qu- 
Dit, i.e. a D-dimensional system. The increase of noise 
level (dilation factor), the LNE and the fidelity are given 
by exactly the same expressions as in the case of qubits 
(see Eqs. I147ll61j) . where the definition of LNE has to 
be adopted accordingly (local white noise of individual 
qudits rather than qubits). 

Similarly, the simulation of pairwise interacting D- 
level systems requires m-body interactions of qudits, 
where 



m 



; log_D 
J \ogd 



(165) 



The simulation of many-body interactions of quDits, say 
n-body interactions, requires effective m-body interac- 
tions of qudits with m = \n log D j log d~\ . 



IX. UNITARY CONJUGATION AND GRAPH 
STATE ENCODING 

We now turn to the second method to generate many- 
body Hamiltonians, unitary conjugation and graph state 
encoding (see Sec. IIIip . In this case, certain two-qubit 
gates U, W (corresponding to a graph state encoding) 
are applied before and after the evolution with respect 
to a single- or two-body Hamiltonian, which leads to an 
effective evolution with respect to a many-body Hamil- 
tonian. The choice of the encoding determines the many- 
body interaction that is generated. In this method, errors 
occur due to 

(i) imperfect application of two qubit gates U, W 

(ii) imperfect application of single and two-body 
Hamiltonians. 

While errors do not appear in higher order 8t' as in the 
commutator method, they can still be significantly am- 
plified. In particular, any simulation involves application 
of gates U and W and hence any error associated with 
these gates. As these gates require the application of a 
basic Hamiltonian H = a z ® a z for time t — 7r/4, errors 
will be of order exp(— kqtt/A), even if 8t' <gC 1. 



A. Fixed graph state encoding 

We will first consider a fixed graph state encoding, i.e. 
the operations U, W only have to be applied at the be- 
ginning and end of simulation process. We start by an- 
alytically estimating the effect of noise in this case, and 
then consider exact numerical simulations which confirm 
our estimates. 



1. Three-body interaction 

The easiest way to generate a three-body interaction 
using a graph state encoding is to consider a three qubit 
system and an encoding operation 



(166) 



where Upg — diag(l, 1, 1, — 1) is a phase gate. With 
this encoding, single qubit operations of the form 
exp(— iSt'a x ) allow one to generate an evolution with 
respect to the three-body Hamiltonian H = ®<J^ <8> 
ai 3 ^ for time t = 6t', i.e. 



U t (St') = Ue 



[/t — g-wStV^Wf'glcrf > (167) 



Phase noise: We now turn to a simple example to 
illustrate how noise effects this process. We consider 
a basic evolution that is governed by the Hamiltonian 
H = <j z ® a z , and where the noise part is described by 
a Liouvillian that corresponds to single qubit dephasing, 
i.e. is given by Eq. 11481 Since noise part and Hamil- 
tonian part commute, the solution of the corresponding 
master equation is given by 



p(t) = V^( P )V^{p)UpU\ 



(168) 



o z 2 ^) and T>^ is a dephas- 



where U = cxp (—ito~ z 
ing map with parameter p — p(t) — exp(— 7ot) (see Eq. 
I104[) . In the remainder of this section we will drop the 
t subscript from the dephasing and depolarizing maps, 

so, unless otherwise stated, = V^ k \ Applying this 
noisy evolution for time t = 7r/4, together with local ro- 
tations among the z-axis (i.e. operations of the form 
exp(— 5ta z ), a noisy phase gate can be generated 

E PGP = £>« {p)V^ (p)U PGP U PG , (169) 

where p — p(t) — exp(— 7o7r/4). The the total (noisy) 
process to generate the evolution with respect to the 
three-body Hamiltonian H for time t = 5t' is thus given 
by 

e P = e$§e™WW)e$3e$8p, (iro) 

where uP(St')p = U x (8t')pU x {5t')^ with U x (5t') = 
exp(—i5t'a x ). Using that the noise process commutes 
with a z , and that V(p)V{p) = V(p 2 ), one can simplify 
this expression and obtains 

£p = V^(p 2 )U t (St')V^(p 2 )V^(p 2 )V^(p 2 )p, (171) 

where Ut(St') is given in Eq. 11671 Notice that T>( x > and 
Ut(St') do not commute. Using Eq. I171l and the chaining 
inequality (Eq. [87| we can derive an exact bound on 
the Jamiolkowski distance and the corresponding fidelity. 
The distance between £ and Ut(St') can be expressed as, 

D(£,U t (6t>)) 

= D{V (1 \p 2 )Ut(St')V (1 \p 2 )V (2) {p 2 )V {3) (p 2 ),U t (St')) 

(172) 
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Applying the chaining property of the Jamiolkowski dis- 
tance we find, 

D(£,U t (St')) 

<D(vW(p 2 )U t (6t'),U t (6t')) 

+D(V^(p 2 )V (2 \p 2 )V^(p 2 ),l) 

= (p 2 ), 1) + D(V^ (p 2 )V^ (p 2 )V^ (p 2 ), 1). 

(173) 

Noting the fidelity of a one qubit dephasing channel with 
the identity is easily calculated and one finds that, 

F^ 1 V), 1) = (174) 

and 

F{V^{p 2 )rt 2 Xp 2 )V^\p 2 ll)= {^AP) ■ (175) 
Which yields the corresponding distances, 

D{V^{p 2 ),l) = \Jl- l -±f, (176) 



D{V^\p 2 )V^{p 2 )V^{p 2 \t) = yi - (I±^(177) 

Thus, D(£, U t (St')) has the upper bound, 
D(£,U t (dt')) 

(178) 

which provides an lower bound on the fidelity, 
D(£,U t (6t')) by 

F{£, Ut(St')) >1- 




(179) 

For small St' a much simpler expression for the fidelity 
can be found that is tighter than the above bound. Using 
that V{p)U x {8t') = U x (St')V(p), where 

T>{p)p = pp+^-(p + ApA^, (180) 

and A = U t (St') 1i a z U t (St'). It follows that the total pro- 
cess can also be written as, 

£ P = ^ t ((^ (1 V)^ (1 VP (2 V)^ (3 VV, (181) 



i.e. a sequence of certain local noise processes followed 
by the ideal evolution. The Jamiolkowski fidelity of the 
total map £ with respect to the desired evolution U t (St'), 
F(£ ,Ut(St')) is thus the same as the one of the map 
£p = V i - 1 "i{p 2 )V'~ 1 \p 2 )V { - 2 \p 2 )V^{p 2 )p with respect to 
the identity, F(£, 1). For small St' , one has that P^P, 
and hence 

£ p » (p^V^ (p 2 )V^ (p 2 )p. (182) 

The fidelity F(£ , 1), and hence F(£ , U t {St')) can be easily 
calculated, and one obtains 

« 1-7T7Q, (183) 

where p = exp(— 7 7r/4) and the last simplification only 
holds for 7o7r/4 <C 1. Under the same approximation the 
exact lower bound on the fidelity that we found in Eq. 
11791 can be approximated to 



F{£,U t {St'))>\- 




(184) 



demonstrating that our bound provides a relatively good 
approximation to the achievable fidelity in this parameter 
regime. 

The local noise equivalent (see Sec. I VI Al for definition) 
can be evaluated under the previous assumptions. We 
find that F = ( i+oxpl-M^ y ~ i _ | 5t / 7] and hence 
the local noise equivalent 7 is given by 

7-^70. (185) 

White noise: Similar arguments can be used to esti- 
mate the effect of white noise, i.e. when the Liouvillian 
describing the noise process in the basic two-qubit in- 
teraction is a sum of local terms described by Eq. 11031 
with s = 1/2 and B — C = kq. The solution of this 
master equation (when taking interaction into account) 
does in general not lead a simple form such as Eq. 11861 
where the ideal evolution is simply composed by noise 
channels. However, one may still use such a form, and 
still obtains a good approximation, in particular if kq is 
relatively small. 

In particular, we have 

p(t) « M (1) (p)M {2) (p)U P U\ (186) 

where U — exp(— ita^ ® o~^) and M is a depolariz- 
ing map with parameter p = p{t) = exp(— not) (see Eq. 
I105|) . As the depolarizing channel does not commute 
with U the best analytic bound that we can find on the 
Jamiolkowski distance is that given by the simple bound 
derived in Eq. [55] Noting that, 

F(M {1) (p)M (2) (p)U,U) = F(M {1 \p)M {2) ( P ),l) 




(187) 
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and that in order to simulate Ut(St') we must use U four 
times we find, 



D{£ , U t (St')) < Wl 



l + 3p 



(188) 



Thus, the corresponding Jamiolkowski fidelity is lower 
bounded by: 



F(£,U t (St')) > 1-16 1 



l + 3p N 



(189) 



Following the same line of argument as for the dephas- 
ing calculation above, we can find a simple approximation 
that beats this bound when St' is small. One finds that 
the total (noisy) process is approximately described by 



£ P » M^(p A )M^(p 2 )M {3 \p 2 )p. 



(190) 



This leads to a final fidelity of the total process given by 



F(£,U t (6t')) 



l + 3p i \ fl + 3p 

T 

3tt 



-«o 



(191) 



where here p = exp(— kqtt/A) and the last simplification 
only holds for Ko7r/4 < 1. If we alternately used the 
exact error bound of Eq. 11891 in the same parameter 
regime we find 



F(£,U t (St')) > 1-37TKO, 



(192) 



which is close to the approximation of Eq. 11911 

Again, the local noise equivalent can be evaluated, 
where in the case of depolarizing noise one can estimate 



F 



4-3 cxp( — n5t ) 



i 



noise equivalent 



2tt 

m 



7«o- 



which leads to a local 



(193) 



Notice in particular that the final expression for the fi- 
delity is independent of 5t' , the time for which the three- 
body interaction should be applied. Such a constant fi- 
delity in turn implies that the local noise equivalent grows 
with decreasing St' . This can easily be understood by re- 
calling that the reliability parameter p for dephasing or 
depolarizing noise decreases with time, and hence a con- 
stant fidelity actually means a larger noise level for short 
time evolutions. 

These observations are in contrast to the commuta- 
tor method, where the fidelity F essentially decreases for 
larger St' . A direct comparison shows that for small val- 
ues of St', the commutator method leads much better fi- 
delities. However, when considering evolutions for longer 
times St', e.g. the generation of gates or the simulation of 
the evolution with respect to a more complicated Hamil- 
tonian generated by Hamiltonian simulation techniques, 



the graph state encoding method performs significantly 
better (see table [J) . In particular when considering time 
evolution over several full cycles (multiples of 2tt), the 
advantage of graph state encoding method becomes ob- 
vious. In this case, the local noise equivalent can even 
be smaller than kq. The latter processes appear, for in- 
stance, when simulating an evolution with respect to a 
time dependent Hamiltonian in such a way that adiabatic 
passage to the ground state of the final Hamiltonian oc- 
curs, and hence the ground state of this Hamiltonian is 
generated [3 [II]. 

We remark that in our estimation of the total fidelity 
we have assumed small values St'. It turns out, how- 
ever, that even when taking all errors into account and 
performing numerical simulations of the total process, 
there is essentially no dependence on St'. In addition, 
our estimates turn out to be very accurate (see table 
|T|. This table shows the resulting fidelities when using 
graph state encoding techniques as described above to 
generate unitary evolutions with respect to Hamiltonian 
(i) Hi = a^P a z 2 ^ a z 3 ^ for time t = 7r/4. In this case, the 
graph state encoding Gi corresponds to U = UpgUpQ 

and single-qubit operations of the form exp(— iSt'ax ) 

are used to generate the desired three-body Hamiltonian. 

Also the usage of an alternative graph state encoding G2 
( 12) 

using only U = U PG together with the application of 

(2) (3) 

two-qubit operators exp(— iSt'a x a z ) to generate the 
desired three-qubit Hamiltonian is shown. 

The latter graph state encoding G2 is also used to gen- 
erate time evolution with respect to the Hamiltonian (ii) 

H 2 = aPaPai^ + £Li ^ for time * = ^/ 4 - In 
this case, Hamiltonian simulation techniques need to be 
applied to generate evolutions with respect to the non- 



r (l)„(2) (3) 



commuting terms af af a\"' and j a x for short 
times St m sequentially. Up to an additional local ba- 
sis change (a Hadamard operation Had on qubit 2), the 
first term is generated by the two-qubit Hamiltonian 
H[ = ax & z 3 \ while the second term is produced from 
the Hamiltonian H' 2 = cri 1) ' a z 2) + a{ 2) + cri 3) . The total 

procedure thus involves the application of [/p^Had^, 
followed by sequential application of H[, H' 2 for short 

times 5t m a total of (n/4)/St m times, and a final appli- 

(2) ( 12) 

cation of Had v U PG . Notice that the same examples 
are considered in the case of the commutator method in 
Sec. IVIIIDI 



2. Many-body interaction 

It is straightforward to obtain methods to generate a 
m-body Hamiltonian using graph state encoding, and 
analytic estimates for the influence of noise. In particu- 
lar, the graph state encoding 



(ik) 

PG ' 



(194) 
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Hi, anal. 


Hi, Gi 


Hi,G 2 


Hi, C 


H2, G '2 


p=0.9 


F=0.8552 


F=0.8545 


F=0.8887 


F=0.7062 


F=0.8545 


p=0.99 


F=0.9851 


F=0.9850 


F=0.9888 


F=0.9305 


F=0.9850 


p=0.999 


F=0.9985 


F=0.9985 


F=0.9989 


F=0.9841 


F=0.9985 



TABLE I: Table shows fidelities to generate gate Ui = 
exp(— in/4Hi) using graph state encoding (column 1: ana- 
lytic estimate using graph state encoding Gi, columns two 
and three: numerical simulation using graph state encoding 
Gi and G2 respectively) or commutator method with optimal 
time Stm (column four) . The last column gives the fidelity to 
generate the time evolution U2 = exp(—iir/4H2) using graph 
state encoding G2 (see text for details). A white noise param- 
eter p = exp(— Ko7r) is assumed that corresponds to a coupling 
strength kq to the thermal bath that would lead to depolar- 
izing map with parameter p when applied for time t — n. 
Notice that this corresponds to p' = p 1 ^ when using the for- 
mula for analytic estimate of fidelity, Eq. 11911 where noisy 
evolution is only applied for time t = n/4. 



together with the evolution with respect to the single- 
qubit Hamiltonian H = for time t = St' allows one 
to generate the m-qubit gate exp(— iSt' ®™ =2 o* ® a - ^), 
i.e. 



Ue^p(-iSt'a^)U^ 



exp[ 



-^'(®? =2< 7W)®4 1 )il95) 



Following precisely the same line of reasoning as in the 
case of three-qubit interaction, m = 3, one obtains in 
the case of phase noise that the total noisy evolution is 



described by (compare with Eq. 1 18 1[) 

m 

£p = U t {5t')V {1 \p m ~ 1 )V^{p m - 1 ) Y[ p( fc V)P-(196) 



fe=2 



We find that the analytic bound on the fidelity of this 
operation given to us by the chaining property of the 
Jamiolkowski distance to be, 

F(£,U t (St')) >1- 



1 +p n 



l+P n 



l+P 2 



(197) 

When St' is small a simple estimation of fidelity can be 
made 



F(£,U t (St')) 



1+p 2(m-l)\ fl +p 



l-2(m-l)- 7o 



(198) 



where p = exp(— 7o7r/4) and the last simplification is 
valid for tt/Ajq <C 1. For the local noise equivalent we 



Similarly, in the case of depolarizing (or white) noise, 
one finds that the fidelity can be lower bounded by 



F(£, U t {St')) > 1 — (to - l) 2 



1 - 



l + ip 



.(199) 



We can also directly evaluate the fidelity in the small St' 
limit to be 



F(£,U t (St')) 



l + 3p 2(m-l)N /I +3 p2x»< < 



1 — 3(m — 1) — Kq 



(200) 



where here p = exp(— Ko7r/4) and the last simplification 
only holds for K()7r/4 *C 1. Again, the final expressions 
are independent of the total time t = St' . The local noise 
equivalent can be determined to be k ~ (m—i)n , 



,61' 



-Kq. 



3. Timing errors 

In this subsection we consider the effect of timing errors 
on the fixed graph state encoding method. Recall that 
111 Subsection En] we introduced the idea of a random 
timing error. Such an error arises due to our inability to 
apply any Hamiltonian for an absolutely precise amount 
of time. 

Consider the 3-qubit simulation protocol of the previ- 
ous subsection as given by Eq. 11671 



U t (St') 



- 7-f( 12 )/-f( 13 ) f ,-«t'4 1 )T f (12)r r (13) 

— u pG u pG e u pG u pG 



Each of the phase gates, Upg, can be generated by 
evolving a two qubit Ising interaction Hamiltonian, H = 
Cz ® Gzi for a time t — 7r/4 and applying single qubit a z 
rotation to each qubit. If we have a timing error on an 
Ising Hamiltonian that generates a phase gate, then we 
can represent this noisy process by the following CP map 



£pg{p) = TUpgpUpg 



(201) 



where T is the CP map given in Eq. 11221 The map T 
is a correlated two-qubit dephasing channel. The noisy 
version this simulation protocol is expressed by 

£{ P ) = r^r^u t {8t')r^T^( P ). (202) 

using the chaining property of the Jamiolkowski distance 
it is possible to bound the fidelity of this operation. Not- 
ing that neither T\2 nor commute with Ut(5t') and 
applying the chaining inequality we find, 



D(£,U t (St')) < 2D(T {12 ^T i - l3 \l). 



(203) 



Assuming that the timing error was Gaussian distributed 
with a standard deviation of a (as was done in Subsection 
IVHl) we find 



F{T^T^\l) = 



1 



(204) 
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where q — e~ . The resulting lower bound on the fi- 
delity is, 



F(£,U t {5t'))>l-A 1- 



(205) 



If we assume that the timing error is very small, that is 
it 2 «1 then we can find a simpler form for this bound. 
Noting that for small a 2 



F(T( 12 )T( 13 ),l)«l-2(7 2 , 
then it follows that 

F{£,U t {8t')) > 1-8ct 2 . 



(206) 



(207) 



The above equations can be simply generalized to the 
m qubit case. Recall that in the m-qubit case that we 
must apply the following encoding operation, 



u e = n 4 ik) 



' PG 



(208) 



k=2 



The result of timing errors on the phase gates in U e is 
given by the following map, 



(209) 



fc=2 



The fidelity, F(£ e ,U e ) = F(U™ =2 T( lk \T), and it 
given by: 



F(£ e ,U e ) = F(l[T( lk \l) = 



k=2 



(210) 



This results in the following error bound for the entire 
simulation protocol, 



F{£,U t {8t')) > 1 - 4(1 - F(£ e , U e )) 



m — 1 ^ 



= I - 4 | I - | l±i ) 1 . (211) 



When a 2 <C 1 the fidelity of performing an encoding 
operation is given by, 

m 

F(£ e , U e ) = F(H T( lfe ), 1) « 1 - (m - l)a 2 (212) 

k=2 

which results in the following bound on the entire simu- 
lation protocol, 



F(£,U t (5t')) > 1 -4(m- l)a 2 



(213) 



So far we have considered the situation where the oper- 
ation UpijUpG is generated by single qubit <j z rotations 
and evolutions of the Hamiltonians i?( 12 ) = and 
_ff( 13 ) = ai^az. We could have in principle generated 



the product of these phase gates by applying single-qubit 
a z rotations and evolving by the collective Hamiltonian 
£f(i23) _ tj(12) _|_ #(13) _ s uc h an evolution would give 
rise to the following CP map, 

£(»3)(p) = T^U P ^U P ^pU^U^. (214) 

Where the CP map T^ 123 ) is a result of timing error on 
the Hamiltonian i/123. The form of 7~( 123 ) was derived 
in Subsection IVIII for the case where the timing error is 
Gaussian distributed with a standard deviation of a (this 
is given in Eq. I130|) . Noting that 



F(£<-™\UMU™) = F(£^ 23 \e-^ H<123) ] 
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the fidelity of this channel was derived and we found it 
to be 



F(T( 123 ),1) 



1 
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Using similar methods to that which we have seen 
above, we can calculate a bound on the fidelity 
F(£, Ut(St' j) given that we used a collective Ising Hamil- 
tonian H ( 123 ) as our entangling resource, 

F(£, U t {St')) > 1 - 4 (l - F(T (123 \ 1)) . (217) 

When er 2 -C 1, we see that 

F(T( 123) ,1) m l-2a 2 (218) 
which results in a lower bound for F(£ , Ut(6t')), 

F{£,U t {5t'))>\-%<7 2 . (219) 

Interestingly, we see that this is the same fidelity that we 
derived in the case where our entangling operations were 
generated by two-qubit Hamiltonian interactions. 

Generalizing to the m-qubit case, we recall that 
we must apply an m-qubit encoding operation, U e = 
Yik=2 Upc ■ This m-qubit encoding operation can be 
performed by single qubit a z operations and an evolu- 
tion for a time 7r/4 of the Hamiltonian, 
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fc=2 



H e is a sum of commuting unitary operations, in Sub- 
section I VII 21 we analyzed the effect of timing errors on 
such Hamiltonians. Because there are only m — 1 
terms in H ei and because of the geometry of the inter- 
action graph, it is possible to find an analytic expression 
for a timing error in the application of the Hamiltonian 
H e . We found the resulting fidelity to be, 



F(£ e ,U e ) = F(T e ,l) 
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where n = m — 1. As we saw above, this allows us t 
bound the fidelity of the entire simulation protocol t 
be: 

F(£, U t (St')) > l-4(l- F(T( 12 - m \lj) . (225 

All of the noise that appears in our analysis of th 
graph state encoding protocol appears within the tin 
ing error of the encoding operation, U e . As we saw, th 
success of the protocol depends on the fidelity F(£ e , U e 
The higher this fidelity, the higher the fidelity of the ove: 
all protocol. Our analysis begs the obvious question; 2 
the presence of timing errors is it better to use a collecth 
Ising Hamiltonian acting on many qubits or many appl 
cations of two-qubit Ising Hamiltonians? We addressed 
more general version of this question in Subsection I VII fl 
Applying the results of Subsection IVII 21 to the case c 
applying the encoding operation U e we find 

m 

F(Y[ T m , 1) < F{T {12 - m \t), (223) 

where we assuming that the standard deviation of the 
timing errors in the two-qubit Hamiltonians is the same 
as that for the collective Hamiltonian H e . Under these 
assumptions using a collective Ising Hamiltonian results 
in a larger fidelity for the simulation. In Figure[8]we have 
plotted the fidelity of performing the encoding operation 
U e using a collective m — qubit Ising Hamiltonian and 
using two-qubit Ising Hamiltonians. We see that for a — > 
that both methods converge to give a fidelity of one. 
As a increases we see that the fidelity of the two methods 
diverges and that this effect becomes more pronounced 
as the number of qubits coupled by U e increases. 

B. Variable graph state encoding 

In our previous consideration, we have assumed that 
a fixed graph state encoding is sufficient to generate 
the desired many-body Hamiltonian using only single- 
and two-qubit interactions. In Sec. IIHI we have dis- 
cussed several examples where this is indeed the case, 
and the strength of the graph state encoding method 
is certainly in such situations. In particular, if the total 
time for which an evolution with respect to a many-body 
Hamiltonian should be generated is large, the graph state 
encoding method performs significantly better than the 
commutator method. 

In certain situations, however, a fixed graph state 
encoding is not sufficient and the encoding has to be 
changed frequently. Any change of encoding requires the 
application of phase gates, i.e. applying the two-body 
interaction described by the Hamiltonain H = o~ z <g> o~ z 
for time t = 7r/4. Such a situation is discussed in detail 
in Sec. IIII E 21 Recall that for a given interaction Hamil- 
tonian, it can be that not all required (non-commuting) 
interaction terms Hi , . . . , H n can be generated using a 
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FIG. 8: (Color online) In these figure we see how the fidelity of 
performing the encoding operation U e changes with the stan- 
dard deviation of a timing error, a. We have examined how 
this fidelity changes with the number of qubits, m, encoded 
by U e . The solid line represents the fidelity of performing 
the encoding operation with a single collective Ising Hamilto- 
nian H e . The dashed line represents the fidelity as a result of 
using a number of two-qubit Ising Hamiltonians Hik to gen- 
erate U e . We see that as a — *■ the two methods give the 
same fidelity. However, it appears for larger a the collective 
method is superior. This effect becomes more apparent as m 
increases. 

fixed graph state encoding. Hence Hamiltonian simula- 
tion techniques need to be applied, requiring the sequen- 
tial evolution with respect to Hamiltonians Hi , . . . , H n 
for short times St' that are generated with help of cer- 
tain graph state encodings Gi, G2, ■ ■ ■ , G n . As discussed 
in detail above, any change of graph state encoding intro- 
duces a certain (fixed) amount of noise. This is reflected 
in a constant fidelity of the simulation process using a 
fixed graph state encoding, or alternatively by an increase 
of the local noise equivalent that is essentially given by 

k/k « n/St'. (224) 

On the one hand, St' needs to be chosen sufficiently 
small to ensure that errors due to taylor expansion in 
the simulation process remain small. On the other hand, 
smaller values of St' lead a larger local noise equivalent 
due to frequent change of graph state encoding. Similar 
as in the commutator method, a compromise between 
these two competing requirements need to be chosen, 
such that for the simulation with respect to a certain 
Hamiltonian for time t the resulting fidelity is maximized. 

X. TELEPORTATION BASED METHODS 

We now turn to teleportation based methods to gener- 
ate many-body Hamiltonians. Recall that the basic idea 
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is to use two-body interactions to generate certain en- 
tangled multi-party states \i(>e)- These states are then 
used as a resource to generate by means of (joint) mea- 
surements the desired evolution with respect to a many- 
body Hamiltonian for time St' . Depending on the chosen 
scheme, the states \i/je) are (i) weakly entangled and se- 
quences of probabilistic teleportation processes together 
with preparation of states with increasing amount of en- 
tanglement are required to make the process determinis- 
tic; or (ii) maximally entangled and involve also an an- 
cilla particle, which is used to chose the evolution time 
St' and to make the process deterministic. 

If two-body interactions are noisy, then the generated 
multi-party states are mixed. The fidelity of the state is 
determined by the quality of the two-body interaction. 
In turn, the fidelity of the multi-party states determines 
the noise level of the many-body interaction. We remark 
that for the approach (i), it is not always clear how to 
use two-body interactions for short times to generate the 
desired state. Even if this is possible, it turns out that ap- 
proach (i) yields smaller fidelities than approach (ii). We 
will thus only consider approach (ii) in the following, and 
assume that maximally entangled states of several par- 
ticles (including ancilla particles) are generated by two- 
body interactions, and processed by measurements. As in 
the case of graph state encoding, here we are again faced 
with a fixed amount of noise, independent of the desired 
interaction time St' . On the one hand, this results from 
the generation of the maximally entangled state \iPe), 
which involves application of two-body interactions for 
times St « 0(tt). On the other hand, the teleportation 
process consists of joint (Bell) measurements, which may 
also require the application of entangling operations that 
further increase the noise level. 

In the following we will discuss a scenario where higher 
dimensional systems are used to simulate qubit systems 
(or more generally lower dimensional systems). Extra di- 
mensions are used to store the particles corresponding to 
the state \4>e) and at the same time the system degrees 
of freedom which should be processed. In this case, the 
measurements in the teleportation process are local as 
they involve only single systems, and one can argue that 
the noise for such local processes is different (and lower) 
than for two-system interactions. Hence the dominat- 
ing part of noise will come from the preparation of the 
entangled state \i/je)- We will elaborate further on this 
scenario, and will also consider the case where additional 
auxiliary dimensions can be used to perform entangle- 
ment purification to reduce the amount of noise, or in 
the case when single system operations are perfect, to 
eliminate noise completely. 



A. Simulation of low dimensional systems using 
high dimensional systems 

In Sec. IVIII El we have argued that low dimensional 
systems of dimension d can also be used to simulate in- 



teracting high dimensional systems of dimension D. In 
this case, many body Hamiltonians are required, even 
to achieve single-system operations on the simulated D- 
dimensional system. Here we consider the opposite ap- 
proach, where high dimensional systems of dimension D 
are used to simulate low dimensional systems of dimen- 
sion d, and we concentrate on d = 2, i.e. qubits. We will, 
however, not embed the qubits into the full Hilbert space 
of the U-dimensional systems, but rather associate each 
D— level system with a single qubit, i.e. we use only two 
of the D dimensions to represent and store the quantum 
information. The remaining levels are used as auxiliary 
levels in such a way that operations on the the qubit sys- 
tems can be performed in an simplified way, or to increase 
the fidelity of such operations and interactions. 

To be more precise, we consider the case where the 
dimension D is some power of two, D — 2™. In this case, 
the D-dimensional system can be considered as a system 
consisting of m virtual subsystems (virtual qubits), where 
one of these qubits, the storage qubit A%, is used for 
storage of quantum information, while the remaining m— 
I qubits A2, . . . , A m are auxiliary systems. We remark 
that this simply corresponds to a labelling of the D- 
dimensional Hilbert space which we choose for notational 
convenience and to clearly distinguish between levels that 
are used for storage purposes and auxiliary levels. These 
virtual qubits do not have any true physical relevance. 

Interactions between the auxiliary qubits and the stor- 
age qubit, as well as interactions between the auxiliary 
qubits, correspond to single system operations on the D— 
dimensional system. We also consider two-body interac- 
tions between two such D-dimensional systems labelled 
by A and B. Depending on the type of interaction, such 
a two-body interaction can couple the two storage qubits 
A\ and Bx, yielding a pairwise interaction, or two aux- 
iliary qubits Ak and Bk (or combinations of these pro- 
cesses). We will assume that single system operations, 
which we also call local operations, can be performed 
with a high fidelity. The noise of such local processes 
is described by a noise parameter pi (when considering 
gates), or by coupling strengthes 7;, Hi when considering 
single system interactions (i.e. interactions between aux- 
iliary qubits and the storage qubit) that are described 
by a master equation. Interactions between different D— 
dimensional systems A and B, which we also call non- 
local operations, can be performed with lower fidelity. We 
will denote the noise parameter by po, or the coupling 
strengthes by 70, «o- We remark that these assumptions 
are natural in the sense that often single system opera- 
tions, i.e. interactions between different levels of the same 
system, are easier to perform than controlled interactions 
between two (spatially separated) systems. For exam- 
ple, atoms or ions have many potentially usable internal 
states (levels) which can be easily coupled via Raman- 
or Microwave transitions, while interactions between two 
such atoms or ions involves more complicated processes 
(e.g. induced dipole couplings, sequences of Laser pulses 
and couplings via motional states etc.). Note that a sim- 
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ilar approach was considered in the context of quantum 
computation in Ref . p59| , where the auxiliary levels where 
used to prepare entangled states shared between systems, 
which where then used as resource to implement non- 
local gates. In 59], the exact requirements that such a 
scheme can be applied were formulated, and in the fol- 
lowing we assume that these requirements are met. In 
particular, the schemes we consider requires the realiz- 
ability of the following operations: (i) (non-local) two- 
system interactions that act only on specific virtual sub- 
systems of each particle, e.g., A 2 ,B 2 , without affecting 
other virtual subsystems; (ii) single-system measurement 
on one virtual subsystem without affecting other virtual 
subsystems; (iii) arbitrary unitary operations on one vir- 
tual subsystem; and (iv) (local) interactions between ar- 
bitrary virtual subsystems Aj , Ak . While (i) requires a 
specific type of two-body interactions between two sys- 
tems A and B, requirements (ii)-(iv) are concerned with 
the ability to locally manipulate a single system. 

Similarly as in the context of quantum computation 
(5t|, here we consider the case that auxiliary qubits are 
used to generate and store the entangled states \iPe), and 
teleportation is used to generate time evolutions governed 
by a many-body Hamiltonian for time St' on the stor- 
age qubits. Note that the involved measurements and 
manipulations are local, and only the generation of \iPe) 
requires non-local operations. In addition, if enough aux- 
iliary levels D « 8 — 32 are available, entanglement purifi- 
cation can be applied to increase the fidelity of the noisy 
entangled states generated by noisy two-system interac- 
tions. The states \iPe) are essentially GHZ states (or 
graph states that are locally equivalent to GHZ states), 
and known multiparty entanglement purification proto- 
cols [H, [56[ can be used to purify these states and to 
increase their fidelity. Consequently, many-body interac- 
tion Hamiltonians acting on storage qubits with increased 
fidelity can be generated. 



B. Perfect single— system operations 

In this section we assume that local operations are per- 
fect, i.e. pi = 1 or 7; = ki = 0, and investigate the in- 
fluence of noise resulting from two-system interactions in 
the generation of ^e)- In the case of n-body interactions 
that should be applied on (virtual) qubits A± , B\, . . . , N\ , 
the state |V>e) is a 2n + 1 qubit state given by 

\if> E ) = l/V2(\^ n \0) E + \4>-) 9n \l) E ), (225) 

where the last qubit corresponds to an additional auxil- 
iary system E and the \4> + ) states corresponding to vir- 
tual qubits two and three of the same system, X 2 X%, 
X 6 {A, B, . . . , N}. It follows that \i/je), when consid- 
ered as a n + 1 system state, is up to local, single system 
operations (that include also local CNOT or phase gates 
which we however consider to be noiseless in this section) , 



equivalent to a GHZ-type graph state 

\$ E ) = [l/y/2(\+) 9n \0) E + |-> 8n |l>*)] ® |0)® n . (226) 

The state \iPe) can be generated by applying phase gates 
Upg between the auxiliary particle E and all other par- 
ties, Ox s-fA b N\Up E G '■ The Bell measurements in- 
volved in the teleportation process act on virtual qubits 
one and two of the same system, and are hence local, sin- 
gle system operation (which are considered to be noise- 
less in this section). Hence the only noise in the process 
comes from imperfect preparation of the state \4> E ) which 
involves n two-system gates. 

An alternative scheme that uses a reduced number of 
virtual qubits per system (two rather than three) is the 
following. We use the n + 1 qubit state \iPe) directly, 
and replace the (local) Bell measurements by local phase 
gates that couple the storage qubits A\,B\,...N\ to the 
qubits A 2 , B 2 , ■ ■ ■ N 2 of the state \ip E ), and measurements 
in the cr^-basis on the storage qubits. Together, this again 
leads to transfer of quantum information to the system 
A 2 , . . . , N 2 . An appropriate measurement on the ancilla 
qubit E, where the measurement direction again depends 
on the outcome of all previous a x measurement, finally 
allows one to implement the operation e~ QO * for arbi- 
trary a in the storage particles (now held in A 2 , . . . , N 2 ). 
Additional Pauli operations -depending on the results of 
the measurements- on the final state are required to ad- 
just the local basis. This scheme has the advantage that 
fewer virtual qubits are involved. In fact, this procedure 
is equivalent to the way of obtaining n-qubit phase gates 
in one-way quantum computation put forward by Browne 
and Bricgcl in Ref. 46]. 

We remark that the auxiliary system E may also be 
thought of being an additional virtual qubit of one of the 
systems particles, e.g. ^3 of system A. In this case, fewer 
two-system gates are required, and the influence of noise 
is lower, yielding to higher fidelities. However one of the 
systems needs to have twice as many levels as the other 
systems, and the process becomes non-symmetric. In the 
following, we will not consider this situation, but rather 
the symmetric situation were the auxiliary qubit is given 
by an additional, independent system particle, and we 
use the second scheme involving fewer virtual qubits for 
simplicity. 



1. Three-body interactions 

We start by considering three-body interactions, 
U(St') = e- 5t ' a f 3 , i.e. a = St'. We consider first the 
case of phase noise, and then the case of white noise in 
two-system interactions. 

Phase noise: 

We consider a basic evolution coupling two virtual 
qubits of two different systems that is governed by the 
Hamiltonian H — o~ z <8> er 2 , and where the noise part 
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is described by a Liouvillian that corresponds to sin- 
gle qubit dephasing, i.e. is given by Eq. 11481 Noise 
part and Hamiltonian part commute, and thus the so- 
lution of the corresponding master equation is sim- 
ple. A noisy phase gate can be obtained using addi- 
tional single-qubit z-rotations, and one finds £pcp = 
2}(l)(p)X>(2)(p)[/ PG p[/j, G , where p = pit) = exp(- 70 7r/4) 
and I? is a single-qubit dephasing map specified in Eq. 

una 

Generation of the state \iPe) using such noisy gates 
leads to a mixed state 



PE 



V^ 2 \p)V^ 2 \p)V^ 2 \p)V^ (p 3 )$ E )(i> E I (227) 



Notice that $ E ) = U p A ^ E) U p B ^ E) U p C ^ E)l 



), and 



slightly different noise process on the final system. We 
find 

£ P = M {MB2C2 \p i )M {A2) {p)M {B2) {p)M {C2) {p)p 1 
where 

M {ABC \p) P = pU{8tf)pU{8^ 

with U{8t') = e- &t '°* , and M(p) is a depolarizing map 
specified in Eq. 11051 with p = exp(— kqtt/A). It follows 
that the fidelity of the process with respect to the ideal 
operation U(St') can be estimated to be 



the phase gates commute with the noise processes. Also 
the phase gates C/ p g lA2 \ U p B q B2 \ J7p^ lC2 ^ used to couple 
the state \4> E ) to the storage qubits, and the a x measure- 
ments on storage qubits Ai , Bi , C\ , commute with the 
noise. Noise acting on the ancilla system E transforms 
after the measurement to correlated phase noise acting 
on the remaining qubits A%, B 2 , C%. That is, for an arbi- 
trary input state p, the output state after the sequences 
of measurements and corresponding local correction op- 
erations is given by 

£p = V^ M \p)V^[p)V^ G2 \p)V^ A2B2C ^{p 3 )U{5t')pU{5ty, 
where 
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(230) 



and for the local noise equivalent k one finds for 8t'n /A <C 
1 
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2. Many-body interactions 
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Vabc{p) P = PP+ 1 -^(P + 

and U(8t') = e~ st '°f\ 

The fidelity of the noisy map £ with respect to the 
ideal process U(a) (where Up = UpW), F(£,U(8t')), 
can easily be calculated and one obtains 



F(£,U(St')) 



It is straightforward to perform the analysis for n— 
qubit operations U(St') = e~ St °* . In this case, a n + 1 
qubit GHZ state is generated by noisy two-body interac- 
tions, and processed by performing local operations and 
measurements. In the case of phase noise, one finds a 
fidelity 
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where the last equality only holds for 7o7r/4 <C 1. For 
the local noise equivalent one finds 7 « -^jlo- 

The corresponding expressions for white noise, where 
p = exp(— KQ7r/4), are 



where the last equality only holds for 7o7r/4 <C 1. The 
local noise equivalent 7 can be determined under the pre- 
vious assumptions, and one finds 
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and the local noise equivalent is given by k 



Note that the local noise equivalent is of the same order 
of magnitude as in the case of graph state encoding (see 
Eq. I185p . and the fidelity is constant, independent of St' . 
White noise: 

A similar analysis can be performed for white noise. In 
this case, noise part and unitary part do not commute, 
and the influence of noise on the ancilla qubit leads to a 
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3. Timing errors 

Analyzing timing errors for this teleportation protocol 
is particularly easy and has many similarities to the de- 
phasing case. Like the dephasing case, we will examine 
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the situation where there are no local errors and that 
all noise in the system comes from the generation of the 
GHZ state \iPe), which is generated via phase gates. As 
we saw above, these phase gates can be generated by a 
two qubit Ising interaction Hamiltonian, H = a z ® a z . 
Here we examine the effect of timing errors on this entan- 
gling Hamiltonian. The action of our noisy phase gate 
on an arbitrary state is defined as, 

£ PGP = TU PG pU PG , (234) 

where T is defined in Eq. 11221 The map T is a correlated 
two-qubit dephasing channel. We note that T will com- 
mute with the coherent processes that generate l^), this 
sibilantly simplifies our analysis. For clarity, we will first 
consider the generation of three qubit interaction. If we 
attempt to generate \^Jb) with these gates the resulting 
state will be, 

PE = T^ E ^ B ^T^ E ^\i> E )(i> E \. (235) 

Now, what is the effect of attempting our teleportation- 
based Hamiltonian simulation protocol using p E as a re- 
source? It is easy to see that noise process commutes 
with the a x measurements on the storage qubits A\,B\, 
and C\. The only non-trivial transformation is as a re- 
sult of the projective measurement on qubit E. We saw 
above that, the effect of measuring a E in the basis that 
we have chosen causes a dephasing channel acting on E 
to transform into a three-qubit dephasing channel acting 
on qubits A2, B2, and Ci- If we instead had a correlated 
phase error between E and Ai before applying the neces- 
sary measurement on E we will have applied a phase flip 
on qubit Ai followed by the same three-qubit correlated 
error as above. The effect of this operation on our proto- 
col is to create a two qubit correlated dephasing opera- 
tion on qubits B2 and C2. That is, the timing noise map 
7~(^2£) g e t s transformed under this teleportation proto- 
col into T( B2C2 K Considering this, the operation induced 
by the teleportation protocol on an arbitrary input p is 
given by, 

£p = T^ B ^T { - A2C ^T^ c ^U{St')p. (236) 

Because the noise process T^^T^^T^ 20 ^ com- 
mutes with U(St'), 

d(s, u(st')) = D(r {A2B2 ^r iA2C2 ^r {B2C2 \ i), (237) 

and 

F(S, U(St')) = f(T (A2B2 ^T( A2C2 *>T (B2C2 ^ (238) 

Thus If we now make the assumption that the timing 
error is gaussian distributed, as we did in Section [VIII we 
can calculate the resulting fidelity to be 

W^'))-(^) 3 + (^) 3 . (239) 



Generalizing this analysis for the generation of the n- 
qubit operations like U(St') = e~ l St is straightfor- 
ward. As a resource we consider an n + 1 qubit GHZ 
state that is generated by phase gates as in the previous 
section. Given an arbitrary input p the effect of per- 
forming the teleportation protocol ideally is to output 
the state U(5t')p at qubits A2-B2...-/V2. Performing the 
same analysis as above and noting that of measuring E 
is to transform operations like 7a 2 e into operations like 
Tb 2 c 2 ...n 2 we nn d that the effect of timing errors is to 
induce the following map 

£p = 'J-( A 2B2---N 2 )^-{A 2 C 2 ...N 2 )j-(B 2 C 2 ...N 2 ) ^ust' p.(2A0) 

The resulting fidelity is, 

F(S,U(6t'))=(^y + (^y. (241) 

C. Entanglement purification and the influence of 
local noise 

In a similar way, one can also take the influence of noise 
in local (i.e. single system) operations into account. In 
this case, the coupling of system particles to the auxiliary 
GHZ state is also noisy, leading to a further reduction in 
the fidelity of the final process. The corresponding noise 
processes are described by similar master equations as we 
use for modelling of two-system interactions, however the 
coupling parameters ki or 7; describing the noise level are 
different. As indicated in the discussion at the beginning 
of this section, we will assume that single system noise is 
smaller than noise in two system operations, i.e. ki <C kq 
and 7; <C jo- 
in the case of phase noise, we have shown that the final 
state is given by Eq. 12281 when assuming that local single 
system operations are perfect, while the influence of noisy 
single system operations is essentially covered by ad- 
ditional noise maps maps V^ A2 \pf),T>^(pf ),V^ C2 \pf) 
with pi — exp(— 7;7r/4). Under the assumption that 
li <C 70i the final fidelity and local noise equivalent are 
essentially determined by 70, with small corrections due 
to additional single system noise. A similar situation is 
encountered in the case of depolarizing noise. 

We will not further go into details about the influ- 
ence of local noise in the standard protocol, but rather 
consider the case where multipartite entanglement purifi- 
cation [55l l55l ] is used to produce high-fidelity resource 
states \ipE) and hence to generate many-body interac- 
tions with high fidelity. That is, we assume that systems 
with additional auxiliary levels — which are used for stor- 
age and purification of entangled resource states — are 
available. In this case, noise in local operations deter- 
mines the achievable fidelity of the resource states, and 
hence the final fidelity of the many-body interaction. In 
particular, the resulting fidelity is independent of amount 
of noise in two-system operations and hence independent 
of kq,7o: provided the noise is sufficiently small such 
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that entanglement purification can be successfully ap- 
plied. For perfect single system operations, the tolerable 
amount of noise in the two system operations is up to 
66 % [59] (or slightly smaller when using direct multi- 
particle entanglement purification) when considering de- 
polarizing noise. Non-zero amount of noise in single sys- 
tem operations leads to a slightly smaller error tolerance, 
however as long as pi = exp(— k/7t/4) is sufficiently close 
to one (errors at the order of percent, i.e. pi « 0.99) we 
find that noise in two-system operations of the order of 
several tens of percent can still be tolerated, without in- 
fluencing the final achievable fidelity of the purified state 
and the resulting many-body interaction. 

An exact analytic treatment of the total process is diffi- 
cult due to difficulties in analytically describing multipar- 
ticle entanglement purification protocols and the reach- 
able fidelities of such schemes. We have thus performed 
a numerical simulation taking all errors due to imperfect 
operations into account. The results are presented in the 
following. 



1. Three-body interactions 

For three-body interactions, the required resource 
state is a GHZ state of four qubits, which can be purified 
using the entanglement purification protocol introduced 
in Ref . [55j, |56( . The fidelity of the resulting three-body 
interaction exp(—St'af 3 ) is plotted in Fig. [§]for St' <C 1 
(St' = 7r/2 15 to be precise) and for a three-qubit gate 
with St' = it/A. As can be seen from the plot, the result- 
ing fidelity is almost identical for the two cases, i.e. to a 
large extend independent of St' . The plot shows the de- 
pendence of (one minus) the fidelity on the noise in local 
operation 1 — pi, where pi = exp(— k;7t) here. Notice the 
linear scaling in the log-log plot. 

We remark the local noise equivalent has no useful 
meaning in this context, as the resulting fidelity of the 
many-body interaction is independent of noise in two- 
system operations (the reachable fidelity of entanglement 
purification only depends on noise in local single system 
operations). More precisely, for any noise parameter kq 
(describing coupling strength to environment for two- 
system operations) that is sufficiently small such that 
entanglement can be generated by the resulting two- 
system operation, we obtain the same local noise equiv- 
alent k which only depends on k;, the noise parameter 
for local single system operations. Notice that n can be 
significantly smaller than the initial noise parameter kq. 
For perfect local control operations and sufficiently many 
auxiliary levels, we can even obtain k = 0. 



2. Many-body interactions 

For n-body interactions, the required resource state is 
a GHZ state of n qubits. Fig. [TU] shows the reachable 
fidelity of the n + 1-qubit GHZ states for different values 




FIG. 9: Fidelity of three-body interaction U = exp(-St'af 3 ) 
for time St' = tt/2 15 « 1(T 4 (red, solid) and St' = vr/4 (blue, 
dashed) using teleportation based method. One minus fidelity 
is plotted as a function of noise in local operation 1 — pi , where 
pi — exp(— K;7r) and 1 — pi — corresponds to perfect local 
operations leading to unit fidelity. 



of noise in local control operation 1—pi- For n — 3, 4 also 
the fidelity of resulting n-qubit interaction exp(— St'af n ) 
for St' = tt/2 15 is shown, which is of the same order of 
magnitude as the fidelity of the purified GHZ state. 



XI. COMPARISON OF METHODS 

In this section, we will recall and compare the main 
properties of the different methods to generate many- 
body interactions. We consider many-body interactions 
exp(— iSt'crf n ) generated by the commutator method, 
graph state encoding or teleportation based method with 
entanglement purification. 

In the commutator method, the fidelity increases for 
smaller St' , even though the local noise equivalent in- 
creases. When implementing a gate (i.e. a certain in- 
teraction for time t = 0(1)) by sequentially applying 
interactions exp(— iSt'af ™) generated this way, there is 
an optimal St'. The fidelity of such a gate is significantly 
lower than the fidelity of a many-body interaction ap- 
plied for short time St' -C 1. 

For graph state encoding, in contrast, the fidelity is 
constant and independent of St' . The fidelity is deter- 
mined by applications of noisy two-body interactions for 
times t = 0(1) (the graph state encoding), leading to a 
significantly lower fidelity for small St' as compared to 
the commutator method. The local noise equivalent also 
increases with decreasing St' , faster as for the commuta- 
tor method. However, also gates or evolutions for times 
t 1 can be simulated with same accuracy, leading to 
an advantage of graph state encoding method over com- 
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XII. SUMMARY AND CONCLUSIONS 




FIG. 10: Reachable fidelity to generate n + 1 qubit GHZ 
states using multiparticle entanglement purification (dots), 
and to generate n-body interactions of the form exp(— 5t'af n ) 
from the GHZ states (x) following the protocol outlined in 
main text. One minus fidelity is plotted against n. Data 
points/curves from top to bottom correspond to noise param- 
eter for local operation (1-pj) = 10" 2 , 1CT 3 , 10" 4 , 1(T 5 , 1CT 6 
respectively, where pi — exp(— k;7t). 



mutator method in such situations. 

The teleportation based method making use of aux- 
iliary local degrees of freedom also leads to a constant 
fidelity, almost independent of 5t' , and hence also gates 
with t — 0(1) can be implemented with same fidelity. 
However, for the teleportation based method using en- 
tanglement purification, the fidelity is determined by the 
noise in local single system operations rather than noise 
in two-system interactions, possibly leading to much bet- 
ter accuracies than reachable with graph state encoding 
or the commutator method. Only for 5t' <C 1 the com- 
mutator method may be more accurate, and for St' ^> 1 
(fixed) graph state encoding may be favorable. The ad- 
vantage of the teleportation based method seems to in- 
crease when considering the implementation of n-body 
interactions with increasing n. It follows that the tele- 
portation based method, together with entanglement pu- 
rification, provides a possibility to successfully simulate 
interacting high-dimensional quantum system or systems 
with many-body interactions with high accuracy. This 
can still be seen as an intermediate stage to full scale, uni- 
versal fault tolerant quantum computation, and should 
be significantly easier to implement than the latter. 



We have introduced and investigated several meth- 
ods to generate many-body interactions from two-body 
interactions. We reviewed the standard commutator 
method which is based on usage of higher order terms 
in the Trotter-Suzuki expansion. We have shown how to 
use unitary conjugation or graph state encoding to gener- 
ate large classes of many-body interaction Hamiltonians 
when using a fixed encoding, and arbitrary Hamiltoni- 
ans when allowing for a variable encoding. Several ex- 
amples, including three-body Hamiltonians with phase 
transition, simulation of interacting d-level systems and 
generation of plaquette interactions have been put for- 
ward. Furthermore, we have investigated the usage of 
maximally entangled states produced by two-body in- 
teractions to generate many-body interactions by means 
of teleportation, i.e. using teleportation based methods. 

We have studied the influence of noise, described by 
generic error models, on the simulation process for the 
three schemes discussed above. For long-time simula- 
tions, the fixed graph state encoding method is favorable. 

However, if one has access to auxiliary systems (e.g. 
higher dimensional systems), then one may use them to 
reduce noise. In particular, one can use higher dimen- 
sional systems to simulate lower dimensional systems, 
and use auxiliary levels to perform entanglement purifi- 
cation. This makes the teleportation based method a 
possible solution to overcome, or at least significantly 
reduce, the influence of noise on a simulation process. 
While this method makes use, to a certain extent, of el- 
ements of measurement based quantum computation, it 
still does not make use of the whole blown up machinery 
of a full fault tolerant quantum computer. We believe 
that our considerations are useful with respect to the on- 
going effort to design and eventually use quantum simu- 
lators based on well controllable, but nevertheless noisy 
quantum systems. 
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